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ABSTRACT 

This thesis considers the design and performance of 
adaptive radar receivers in the presence of clutter with slowly 
time varying statistics. 

An optimum radar receiver for detecting a Swerling 1 
target in the presence of clutter is derived. Considering 
the complexity involved in the implementation of such a 
receiver in a changing clutter environment, two suboptimum 
receivers are proposed and expressions for their probability of 
detection performance are derived. These receivers are 
adaptive in that the slowly time varying clutter statistics 
are taken into account through the maximtun entropy method of 
spectral analysis of the received signal. 

Simulation results show tha.t the first of these 
receivers, termed the Struct'ured Optimum Receiver, can operate 
satisfactorily with as few as 16 hits per scan. The performance 
of this receivor which is always superior to that of the 
conventional matched filter receiver is compared with that 
achievable with a received signal without clutter. In the 
wide range of situations considered, the degradation in 
performance is small except for low target Bopplejr shifts 
and large clutter spectral widths. 
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The second receiver, termed the Structured Receiver, 
can easily he implemented relative to the structured optimum 
receiver particularly if noncoherent detection is employed* 
The performance of this receiver under both coherent and 
noncoherent detection modes is compared with that of a 
moving target indicator (MTI) receiver. 



CHAPEER 1 


IRTRQDUCTIOJT 


1 . 1 Back|g:round and Previous Related Work 

In radar systems, the receiver's ability to detect 
target echoes is limited not only hy additive thermal noise 
hut also by another typo of inter fore nee called clutter. 
Clutter is defined as a conglomeration of unwanted radar 
echoes and arises from the hack scattering of electromagnetic 
energy from possibly extended and fluctuating reflectors. 
Clutter signals differ from therml noise in tha.t they are 
caused by the radar's own transmission. Also what constitute 
clutter producing reflectors clearly depends on which class 
of targets a particular radar is intended to detect. In 
general, clutter signals are statistical in nature and the 
statistical modelling of radar returns from ground, sea etc, 
have been the subject of many publications [l-S)* Although 
several models have been proposed, none can be considered to 
adequately explain the observed returns. In a clutter 
model C5| most often used for theoretical analysis, the 
complex envelope of the clutter signal is considered to be 
a sample function of a zero-mean complex Gaussian random 
process. Under wide-sense stationarity and uncorrelated 
scattering ass'umptions , the statistics of clutter may then 
be described by an associated scattering function which 
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gives a measure of the distribution of the returned clutter 
energy in the delay- Doppler domain [53. Although, in 
practice, these statistics change with time, they nia.y be 
modelled as q.ua,si time-invariant. 

Radar targets are also modelled statistically 
and various target models have been proposed for different 
situations by Marcum [6-73 » Swerling [ 8-93 ► V/einstock [103 
and Mitchell [ll). In these models, the radar cross- 
section of a target is treated as a random variable, fhe 
various models differ in their prescription of the 
probability density function and fluctuation rate of the 
radar cross-section of the target. For example, the 
radar cross-section of a Swerling 1 target model is slowly 
fluctuating and has an exponential density function. 

In most of the previous related work on the design 
of radar receivers for clutter rejection, it has been assumed 
that the statistics of clutter are a priori known. Also, in 
addition to receiver design alone, considerable emphasis 
has been placed on the choice of the transmitted waveform 
and joint transmitter/ receiver design. In fact, receiver 
design is very closely linked to signal (or transmitter) 
design. 

There are three principal approaches [5) to the 
design of receivers for processing of radar echoes in the 
combined presence of clutter and thermal noise. The first 
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approach is to find the optimum receiver (optimum in terras 
of a certain specified performance criterion) for a given 
clutter scattering function. However, finding this optimum 
receiver for any arbitrary clutter scattering function is 
extremely difficult. Moreover, the optimum receiver may 
turn out to be too complicated to be realised in practice. 

Another approach is to design a conventional matched 
filter receiver assuming that only white thermal noise is 
present. The degradation in performance due to clutter 
is then minimised by signal design. For any given clutter 
scattering function, this consists of minimising the 
common volume of the scattering function and the shifted 
signal ambiguity function in the delay-Doppler plane. The 
best signal will depend on the target's location in the 
delay~Doppler plane as well as on the scattering function 
of clutter. In general, receivers designed using this 
approach are suboptimum in performance. 

The third approach is to structure the receiver- 
so that it is implementable - and then, within the 
constraints imposed by that structure, optimise the 
performance, of such a receiver. In general, the performance 
of such a receiver will also be suboptimum. However, this 
structured approach is attractive as it permits design of 
a receiver that works better than the conventional matched 
filter receiver (second approach), but can be constrained to 
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bo lesB complicated than the optimum receiver* In addition 
to receiver optimisation using the structured approach/. 

3 oint transmitter/raceiver optimisation can be done for a 
class of transmitter signals/receivers. 

The design of optimum receivers for clutter 
rejection were first inve stilted by Urkowitia [123 and 
Manasse [133* Both assumed a rahge ' invariant scattering 
function with a Doppler spread which was small compared 
to the signal bandwidth* Urkowitz ignored thermal noise 
and proposed the inverse filter for clutter rejection* 

Manasse considered the presence of thermal noise and 
showed the dependence of the signal to interference (thermal 
noise plus clutter) ratio at the output of the receiver 
on the signal bandwidth. He found the best signal to be a 
wideband signal. Vfesterfield et.al* [143 considered the 
approach of designing a conventional matched filter, and 
then optimising its performance in the presence’ of clutter 
by signal design* Stutt and Spafford [153 proposed a 'best* 
mismatched filter or receiver where the signal to clutter rati' 
is maximised. Rummler [16) » Spafford [173 > Delong and 
Hofetetter [18] also investigated the signal and receiver 
design problem for a slowly fluctuating point target 
assuming that the scattering function of clutter is 
a priori known. Most of these investigations were carried 
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out for a claes of transmitter signe.ls /receive re. Delong 
and Hofstetter [18] and Spaffordtl 7 ) alBO-proposed steps-’ for 

r 

j oint optimisation of both the transmitter signal and the 
receiver for a class of transmitter signals/receivers. More 
related work on the design of radar receivers and transmitter 
signals cen be found in [19-253. 

However, in most practical radar systems, because 
of lack of knowledge of the statistics of clutter, the 
moving target indicator (MTI) receiver is used to 
distinguish moving targets in clutter. MTI receivers are 
simple to implement and are based on the reasoning that 
ground and sea clutter returns have a very small Doppler 
spread and can be filtered out or cancelled in delay line 
cancellers. MTI receivers are obviously suboptimum and 
their design ignores the presence of thermal noise, A 
variation of the MTI called the moving target detector 
(MTD) has been developed at the Massachusetts Institute 
of Technology [26-273. Ik. the MTD, fast Fourier transform 
(PFT) techniques are used to estimate the spectrum of the 
received signal and filter out clutter returns. 

In this thesis, an optimum radar receiver for 
detecting a Swerling 1 target in the presence of clutter 
is obtained and it is found that implementing this 
receiver in a changing clutter environment is extremely 
difficult. In view of this, this thesis primarily deals 
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with, the design and sttidy of two adaptive structured radar 
receivers for the detection of a Bwerling 1 target in the 
presence of clutter. These receivers are adaptive in the 
sense that their parameters are adjusted or updated 
in response to changes in the statistics of clutter. The 
quasi time~inva, riant nature of the statistics of clutter 
is taken into account through the maximum entropy method of 
spectral analysis of the received signal during those scans 
when there is no target or a decision to that effect is 
made , 

1*2 Organisation of the Thesis 

This thesis is organised as follows. In this thesis 
the various bandpass signals, systems and random processes 
are described in terms of complex envelope notation [5]. 

In. Chapter 2, the radar target detection problem 
in the presence of clutter is posed as a hypothesis testing 
problem. Various target models in the context of pulsed 
radars are considered. The complex envelope of the clutter 
signal is modelled as a zero- me an complex Gaussian random 
process. Some of the receiver structures obtained in later 
chapters as solutions to the target detection problem posed 
in this chapter, require spectral estimators as subsystems 
in order to be adaptive to the changing clutter environment. 
As the maximum 0 ntropy Jff.ethod (MBM) of spectral estimation 
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is most non-committal in its processing of the ohservations 
made, we have considered it appropriate to use the MEM for 
estimation of the required statistics of clutter, Accordingly, 
in Chapter 2» the MBiM of spectral analysis of a stations.ry 
time Bv;3ries is hriefly reviewed. 

In Chapter 3, the optimum receiver for Sworling 1 
target models and expressions for its performance are 
derived. It is found that optimum receivers which 
continually estimate the required clutter statistics and 
then update their parameters are difficult to implement 
with present day technology. In view of ■ this# suhoptimum 
receivers >which are simpler to implement and whose parameters 
can he determined adaptively using the MEM, are proposed and 
evaluated in Chapters 4 and 5. 

In Chapter 4> a structure is imposed on the 
receiver and within the constraints imposed by that 
structure, an optimum receiver is obtained for Bwerling 
1 target models. Expressions for the performance of such 
a roceiver, termed Structured Optimum Receiver, are 
derived. A study is made of how the MBIM can he used to 
estimate the required statistics of cluttor during those 
scans when no target is present in the range hin of interest, 
thereby enabling the parameters of the structured optimum 
receiver to be determined adaptively, The performance of 
such a receiver is evaluated by computer simulation and 
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numerical computation for both the cases when the required 
clutter statistics are a' priori known and when the MEM is used 
to estimate the required statistics of clutter* These 
performances are compared with the performance of a 
conventional matched filter receiver, A large number of 
receiver operating characteristic (ROC) graphs indicating the 
performance of the struct^^red optimxmi receiver for varying 
conditions of clutter power, clutter spectral width, target 
Doppler shift etc, are including in this chapter. 

In Chapter 5, an adaptive sub optimum re ceiver» termed 
Structured Receiver , which in many instances is simplor to 
implement than the structured optimum receiver, is proposed. 
This receiver is based on interpreting MEM as an attempt to 
determine the structure of a whitening filter for a time 
series. This chapter also discusses how the structured 
receiver can be used for both coherent and noncoherent 
detection. Expressions for the performance of structured 
receigers for both coherent and non-coheront detection modes 
are derived. The performance of these receivers are evalua.tod 
by computer simulation and nimierical computation and compared 
with the performance of a double delay line canceller MTI 
receiver. Again a large number of ROC graphs indicating the 
performance of these structured receivers for varying 
conditions of clutter power# target Doppler shift, 
clutter spectral width etc. are included in this chapter. 



Chapter 6 concludes the thesis with some suggestions 
for future work. 

Appendix A describes the computer simulation scheme 
used in the performance evaluation of the structured optimum 
and structured receivers. Appendix B gives a listing of 
the main program and the subroutines used in the performance 
evaluation of the structured optimum, structured, matched 
filter and MTI receivers. 
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CHAPOJiR 2 

THE) RAD^m MiiljjiOlIOH IROBiaiM AM) A RliiVIBW OF TWi 
MAXIMUM EFTROPY FISiTHOD OF SPECTRAL AMYSIS 

The first part of this chapter formulates the radar 
target detection problem in the presence of clutter and 
thermal noise. Section 2,1 considers the radar detection 
problem as a binary hypothesis testing problem. Section 2*2 
enumerates various target models in the context of pulsed 
radars. This is followed by a discussion on the modelling 
of clutter in Section 2*3* Tii® pulsed radar detection problem 
is stAted in mathematical terms in Section 2*4. 

Some of the receiver structures obtained in later 
chapters as solutions to the radar detection problem posed 
here, require spectral estimators as subsystems* As 
maximum entropy method (MiilM) of spectral analysis is most 
non-committal in its processing of the observations made, 
we have considered it appropriate to use the MEM for 
estimation of clutter spectra. Accordingly, in the second 
part of this chapter, we briefly review the maximum 
entropy method of spectral analysis of a stationary time 
series. 

2*1 Detection as a Hypothesis Testing*^.- Problem 

A radar* target detection problem may be posed as 
a binary hypothesis testing problem, which may be stated 
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as one of finding a decision function of the observations 


E(t) = 

%(t) 



E(t) = 


fjS t < Tj 

= Ho 


which optimises a specified criterion of performance. 
contains both signal and noise while RQ(t) contains only 
noise. The nature of assumptions made concerning the 
statistical information of the received waveforms on 
hypotheses and determines the complexity of the 
decision function in respect of cost, ease of implementation, 
sensitivity to small changes in the assumed apriori 
information and performance. In radar it is difficult to 
assign a priori probabilities for the occurrence of or 
Hq, nor is it possible to assign realistic cost functions 
for the decisions made. Therefore the Neyman-Pearson 
criterion is most often used in arriving at a decision 
function* This criterion maximises the probability of 
detection of a target (P^) subject to the constraint that 
the probability of false alarm (Rp) is equal to or less than 
a certain prespecified number. It can be shown [28-30] that 
this criterion leads to the likelihood ratio test as the 
optimum test. This test is defined by 




Hi 

- A(R) > X 


( 2 * 1 ) 


where Pj^/h^CR/Ht ) and p^/^ (R/H^) 


R/Hq^ 


are the multidimensional 
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conditional probability density functions of the observations 
given that and are true 'respectively, A (R) ia the 
likelihood ratio and X is the threshold for comparison*. 

Thus the optimum receiver processes the observations 

to find the likelihood ratio A (R) and then compares A(R) 

with a threshold in order to make a decision. If we denote 

the probability density function of A when is true as 

Pa / jj (A/Hq), and when is true as P/v/jj (A/%)» then 
0 1 

00 

* I Pa/h and Pjj * J (A/H^)dA (2*2) 

The threshold X is obtained by setting Pj, equal to the 
maxiraum permissible value of false alarm profcability and 
solving ( 2 . 2 ). 

2 . 2 Target Models for Pulsed Radars 

FumorouD... target models.^ have bo en, proposed for 
different situations [ 6-11}. Marcum [5 ] assumed that 
the target can be modelled as a nonfluctuating or constant 
reflector. However^ in general, the radar cross-section of 
a target is a statistical variable. Swerling in his classic 
paper [ 8 } proposed four different models of target 
fluctuation. These four fluctuation models are commonly 
referred to as Swerling 1, 2 » 3 and 4 models respectively. 
They are as follows J 
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Sweyling; X model 

The amplitude of the entire reflected pulse sequence 
from the target during a particular scan is a single random 
variable with a Rayleigh density function. In addition, 
there will be phase coherence between the various pulses in 
the returned sequence except that the initial phase of the 
sequence will be a statistically independent random variable 
uniformly distributed over 2% radians. The density function 
for the target cross-section (a) will therefore be given by 

p(o) = i exp(-cT/ 0 ) (a > 0) (2.3) 

a 

where a is the average cross-section. In complex envelope 
notation, the complex envolopo of the entire pulse sequence 
is multiplied by a single -zero-mean complex Gaussia.n random 
variable as it gets reflected from the target. However, the 
complex mtiltipliers in different scans are statistically 
independent. This type of fluctuation is called scan to 
scan fluctuation. 

S werllng 2 model 

Here the amplitude of each reflected pulse in the 
sequence is a different statistically independent Rayleigh, 
distributed random variable. In addition, the initial 
phase of each pulse will be a different statistically 
independent uniformly distributed random variable. The 
density function for the target cross-section is again 
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given by (2.3). In complex envelope notation, the complex 
envelope of each pulse in the sequence ie multiplied by 
different statistically independent zero-mean complex 
Gaussian random variables as the sequence gets reflected 
from the target. This type of fluctuation is called pulse 
to pulse fluctuation. 

S werlinp: 3 model 

This is scan to scan fluctuation and the density 
function for the target cross-section (o) is given by 

p(cr) - exp(-2o/ 0 -) (o > 0) (2.4) 

a ^ 

S werling 4 mode l 

This is pulse to pulse fluctuation and the density 
function for the target cross-section is given by (2.4). 

It would be well to indicate in which situations the 
various models may apply. The density function given in (2.3) 
is observed when the target consists of many independent 
scattering elements, of which no single one predominates. 

Most available observational data on aircraft targets 
indicate agreement with the. exponential density. The 
density function given in (2.4) corresponds to that of 
a target having one main scattering element that predominates, 
together with many smaller independent elements. One such 
target may be a spherical satellite with extensions. The 
scan to scan fluctuation model would apply to targets such 
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as jet aircrafts or miss lies, or for radars having a fairly 
high pulse repetition rate, The pulse to pulse fluctuation 
model would apply to propeller driven aircrafts if the 
propellers contribute a 3.arge portion of the echoing area, or 
for radars having a low pul- so repetition rate. 

The above four fluctuation models of Swerling have 
been foimd to bracket a wide range of practical cases, 
particularly for aircraft targets. In 1970, Swerling [93 
extended his original' models to cover more cases of interest 
by considering three families of radar cross-section models; 
the chi-square family, the Rice family and the log-normal 
family. 

In this thesis., however, only Swerling 1 target models 
are considered. Also the targets are assumed to be small 
compared to the radar resolution cell, that is, they are 
point targets. 

2 • 3 Modelling of Clutter Returns 

¥ 

In this thesis, it is assumed that ground, sea and 
weather returns constitute clutter. The statistical modelling 
of ground and sea returns have been the subject of many 
publioations [1-5, 31 , 32,33 3 and although several models 
exist none can be considered to adequately explain observed 
returns. The clutter model [53 most often used for theoretical 
analysis is that of a large number of independent scatterers, 
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with, no single one predominating, within the radar resolution 
cell. The phases of the retxirns from each ecatterer are 
asstuned to he random and uniformly distributed over 2n radians. 
Also each phase angle is statistically independent of the 
other. In general, clutter producing targets or clutter 
targets are extended fluctuating targets. More specifically, 
let f(t) he the complex envelope of the transmitted signal. 
The reflected signal from any range interval ( X, X+ dX ) will 
be the superposition of a large number of reflections with 
random phases. If the phases are statistically independent 
and none of the scatterers predominate, central limit theorem 
based arguments can be used to characterise the reflection 
coefficient as a zero-mean complex Gaussian random variable. 

If the orientation and composition of the reflectors change 
with time, the reflection coefficient must be modelled as 
a random process, The reflected signal from any range interval 
(X , X + dX ) can be written as 

rio(t, X) « f(t -X ) b(t - I, X ) dX (2.5) 

where b(t, X ) is a zero-mean complex Gaussian process whose 
independent variables are both time and space. The return 
from the entire extended target is a superposition of the 
returns from the incremental intervals and is given "by 

n^(t) « f f(t - X ) b(t - I , X ) d X (2.6) 

The complex envelope of the clutter return n^j(t) is a sample 
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function of a zero-mean complex Geuesian process* It can be 
seen that the covariance function of n.(t) is signal dependent. 
Under the assumptions of uncorrelated scattering from different 
range intervals and wide sense etationarity of 'B(t, X)# the 
covariance function of h(t, X ) can be written as 

B[b(t, X ) b* (u, X^)] - Kp^(t-u, X ) 6(t - Xj^) (2.7) 
where 6(t) denotes a imit impulse function and 
B[ ,3 denotes an ensemble average. In actual practice 

l«W 

b(t, X ) will not be wide sense stationary. However* assuming 
b(t| X ) to be wide sense stationary over a short interval of 
time or over several scans is realistic enough. 


Assuming that the clutter extends from zero range to 
range LT seconds* (2.6) can be written as 




f(-t - X) ^(-1 - X ) ax 


( 2 . 8 ) 


Now assume that the complex envelope of the transmitted signal 
is bandlimited to Ng/2 around its carrier frequency* or 

?(f) 0 for |f I > ¥g/2 (2*9) 

where P(f) is the Fourier transform of f(t)» 

1 

Also let T « ^ where W i Wg* The criterion for 
choosing T will be explained later. 

Following Kailath [ 54 3 and Van Trees [53 tiie 
clutter return h^(t) in (2»8) can be written as the output 
of a tapped-delay line filter with taps spaced T seconds 
apart, when the input to the filter is ?(t). Or 
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n (t) « f{t-kT) (2«10) 

fcsO 

The tap gains bjj.(t)’e are 2;erc>-ffiean complex Gaussian random 

*•> 

processes and if X ) in (2*7) is essentially constant 

■with respect to X over an interval T seconds long, then 


E[hi^(t) hJCu)] 


T kT) , k X I 

0 , k ^ A 


( 2 . 11 ) 


In the case of moving weather clutter, the clutter return 
will be centered around a mean Doppler frequency. If to^g, is 
the mean Doppler frequency shift of the returned signal due 
to the motion of weather, then using the above arguments > the 
complex envelope of the clutter signal is given by 



exp(jw^gp t) f(t-kT) 


( 2 . 12 ) 


In this thesis (2.10) will be used to represent the complex 
envelope of the returned clutter signal. However, (2.12) 
could also have been used to represent the complex envelope 
of the returned clutter signal. 

2.4 Mathematical Statement of the Problepi 

let the transmitted signal be given by ^ 

S^(t) = Ifa (2-13) 

oe> 

where w is the carrier frequency and j lf(t) dt = 1. 

^ -*00 

is the energy in the signal, .When a target is at range 
R(t) and moving inwards with velocity V (R(t) » R^-Vt), the 
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received signal from the target is given ty 


S^(t) = V"2 Re[VS)^ B f(t ~ 't(t)) exp[3w^(t-T(t))) 3 (2,14) 

■V7here T:(t) is the round trip delay time, B is a zero-mean 
complex Gaussian random variable if a Swerling 1 or Bwerlang 2 
target is assumed. A signal received at time t would have 
been reflected from the target at time (t~ ('T(t)/2)). At that 
time the target range would be 

R(t - = Rq ~ V(t - •^) (2.15) 


By definition 

2E( t - ^5^) 

•u(t) = ? 

c 

where c is the velocity of light,- 


(2.16) 


Substituting (2,16) in (2.15) and solving for T(t) 
we obtain 


t(t) = 


2Ro/c 
1-f V/c 


2(V/c)t 
1+ V/c 


(2.17) 


Por target velocities of interest > V/c << 1. Then 


T(t) - 2(Ro/c) ~ 2(V/c)t (2.18) 

let - 2 Rq/c be the nominal range delay. The received 
signal is given by 

S^(t) « V2 Re[pl^ B f(t - + ^t)exp[jmQ(t-f ^t-T^)33 

(2.19) 
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If ¥gT<< ■where Wg is the bandwidth of f(t) and is 

the duration of f'(t), we can ignore the compression or 
stretching of the time scale of the complex envelope. The 
Doppler frequency shift is given "by 

Wp = w^(2V/c) (2.20) 

Absorbing the exp (-3 phase term in the phase of B# the 

received signal is given by 

Sp(t) = V2 He[pl^ B f(t - T^) exp(jwQt+j(Opt)3 (2.21) 

The complex envelope of the received signal^ Sj,( t )) from a 
target at nominal range delay is given by 

S;^(t) = B f(t - T^) exp(jW2jt) (2.22) 

If Iil[ I B 1^3 = t the expected value of the received 

signal energy is equal to 2E^a^. The additive receiver 
thermal noise w(t) will be modelled as a zero-mean bandpass 
white Gaussian process['5 3, It can be represented as 

w(t) K f2 Re[w(t) exp(3WQ't)} (2.25) 

where is the carrier frequency and w(t) is the complex 
envelope and is a zero-mean complex Gaussian process. If the 
bandwidth of w(t) is larger than the other bandwidths in the 
system of interest, 'w(t) can be represented as a complex 
white noise process. 

let r(t) be the complex envelope of the received signal 
in the presence of clutter and thormal noise. 
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Ueing the notation developed in the earlier sections^ the 
radar detection problem - to detect a target at range delay 
and Doppler in the presence of clutter and thermal 
noise - can now be stated as the following binary 
hypothesis testing problem: 


r(t) = §^(t) + n^(t) + w(t) 
r(t) ~ + w(t) 



Til- t< 


(2-24) 


Substituting for S^(t) from (2.22) and for n^Ct) from (2.10) 
we get, 

r(t) « B f(t-T^) exp(jti^t) + f b 


f(t) = 


^ )Lo 

+ w(t) ; 

D ^ 

Vil+ l b^(t)l(t-kT) 
+ v(t) 


Ti<-t< If 


: H, 


Also 


= B[bjj.(t)? = E[w(t)3 = 0, B[w(t)w*(u)3 = 11^6(t-u) 

.2 
B 


E[|b|® } == 2e» , B[b|^( t)bl (u)) =:pr Kysrj(t-u, kT), k * a 


0 , k ^ A 

«**» < t,V.-< «*> 

(2.25) 

This hypothesis testing problem will now be modified in the 
f oil owing way : 
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In a particular scan let M- f ( t)» wave forms or pulses 

te transmitted with, a repetition rate of l/T^* Let 

where is the dur^'tion of f(t). We assume that the bj^(t)’s 

are so narrow "band that during the time a pulse - f(t) - is 

illuminating it, the bj^(t)'s are random variables. This means the 

sample functions of hj^(t)'s are so slowly varying that during 

the time a pulse is being reflected ( seconds ), b^(t) 

However, when M pulses are being transmitted during a parti~ 

cular scan, for each pulse, the b^(t)‘s may be different random 

variables. For example, for the first ^(t) , t) « 

the second pulse in the pulse train f(t-Tp), tjj.(t) ’ “ 

so on. It essentially means that the reflection coefficient 
#*>* 

b(t- ^ ) iu any given range interval ( X ,X + dX) will be a 
random variable for the duration T^ when the waveform f(t) is 
being reflected from that interval. In-, the context of pulsed 
radars, and ground, sea and weather clutter, this is a very 
realistic assumption, 

If the complex envelope of the transmitted signal is 
given by 

M-1 

St(-t) = Vlit I - ml ) (2.26) 

m=o ^ 

where 

oo 

J f(t-nTp) f*(t-mTp)dt = 6^ as T^ » 
the hypothesis testing problem becomes, 



Hi 
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M~1 

r(t) = I B i(t-Tr^-mT )exp(jovt) 

rrt=o ^ 


+ 

M-1 

I 

m=: 0 

f 

fc=o 


f ( t-kQUmT ) +w( t ) 

Jr 

'Si 

VEt 

O H 

L 

1&=0 


f( t“'kT-mTp)+w(t) 

'So 


Ti < t <• 


(2.27) 


where are jointly complex Gaussian random variables with 

JStbi,,!,] = 0, E[ 1 

where 6^^ is a Kronecker delta function. (2.28) 

This hypothesis testing problem will again be modified in the 
following way; Let the entire range duration of interest be 
divided into range bine and let the length of each range bin 
be T seconds. Let there be in all (N+l) range bins* Now 
instead of testing for a target at: a particular range delay 
we will test for a target in a particular range bin by 
testing at the centre of the bin. Now T should be small 
enough so tha.t there is little degradation in performance 
even if the target is at ±1!/2 from the centre of a range bin* 

If we let be in the nth range bin, we can replace by nT 
in (2.28). Without loss of generality we will aseume n * o 
and test for a target in the zeroHh range bin, 

¥e also impose the following conditions on the waveform: 


•00 


J f(t-kT~mI^)f*(t--jT-nTp)dt ^ 0 forj,k« 0,..,,N 

, when m ^ n 
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and 


V*/ .v 

J f(t~kT-.mT^) exp(jw^t) f *( t~j CP-nTp )dt = 0 

for j ,k = 0,1, . .. , 1 T 

when m n (2. 29} 

Then the hypothesis testing problem becomes 


r(t) ~ p,^ I B f(t-mT ) exp(j(o t) 
m= 0 .y " 


M^l L 


^ A \k ) + w(t) :] 

m= 0 Ifc 0 ^ 



In Chapters 5 “ 5» various receiver structures are 
obtained as solutions to this hypothesis testing problem. 

Some of them involve the use of the maximum entropy method(r4EM) 



26 


the entropy is defined as 

oo CO 

% = ” J **• ! pCx(i),... x(n)) X 

...CO mmOO 

logp[X(l), ... X(N)3 dX(l), ... dX(U) (2.55) 

where p[ r®P 2 ^® sents the N dimensional joint prohahility 
density of the random variables X(l), ... X(N). Jaynes contends^ 
that the probability density function p[, 3 that is in agreement : 
with the 'testable* available data but is maximally non- committal 
about unavailable information is the one that maximises entropy. ! 
To apply this method to power spectrum estimation, we observe j 
that power spectrum representation implies second order 
characterisation of a stationary random process. Shannon [453 
showed that of all the processes over (- oo, oo) having the same 
first and second order moments, the Gaussian process has the 
maximum entropy. Thus the MEM constrains the process to be i 

modelled as Gaussian. How let Xjq. = [X(l), X(2)>»..» X(H)5^ | 

where the superscript T denotes transpose of a vector or a I 

matrix. Then assuming the process to be real and zero-mean, | 
the joint probability density of these observations can be 
written as 

P(X]j) ~ ' " ' ©xp(— Xjj Tjjj- Xjj) ( 2 . 54 ) 

where D 0 t(,) denotes the determinant of a matrix and 
Tjj- denotes the (H x H) auto covariance matrix given below, 
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R(0) R(l) R(N-.l) 

R(l) R(0) ... R(H*2) 


R(R-l) R(K-2) ... R(0) 

and 


R(p) = B[X(m) X(m-p)3 


(2.35) 


Substituting (2.34) in (2*33) the entropy can be shown 
to be [4 6) 

Hjj * 4- log[Det(Tjj)3 + Constant (2.56) 


For increasing R, diverges# so instead entropy density h 
as defined below is maximised 


h = lim 



(2.37) 


If initially it is assumed that the autocovariance function 
R(k) for I k| < R are known, the MEM strategy results in the 
extrapolation of the autocovariance function R(k) for |k | ^ R 
such that 


= 0 for [kl >, N (2.38) 

8R(k) 

subject to the constraints inposed by the known R(k)'s 
for 1 k| < R, That is, R(k)'s for | k | are chosen so as 
maximize the entropy density h in (2.37) without perturbing 
the autocovariance values already known. The autocovarian.ee 
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matrix Tjj is eymmetric and the elements along any diagonal 
are identical. This Toeplitz structure is'. instrumental in 
establishing a relation bet'ween the power spectral density 
S(f) of the process and the entropy density h. The following 
relation between the determinant of Tj^ and the power spectral 
density S(f) of the process can be established [463» 

lim [Det(Tjj) 3 ^/^ = 2 fj^ exp[{L/ 2 %> f logS(f)df 3 

provided Det(Tjj-) 7 ^ 0 (2*39) 


where fj^ = l/2At and At is the time interval between 
successive observations. After certain manipulations using 
( 2 . 38 ) » (2.37)j (2.38) and (2.39). it has been shown by [463 
that the power spectral density S(f) that maximises the 
entropy density h is given by 


S(f) 


^N-1 


W-l 




1+ 1 exp(-o 2Tcfk A t)l * 


k=:l 


(2.40) 


where the coefficients ^ 2 * “*' '^N -1 %-l obtained 

by solving 

N-1 

R(j) + \ R(d-k) =. 6^0 ; 0 = 0,1...., N-1 

( 2 . 41 ) 

Equation (2.41) can also be written as 
■ N-1 

R(3"k) s ; 4 =* 0,1.. ...N-l (2.42) 

with ss 1, 
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Thus to obtain the power spectral density, only and 

^N-l determined. The Toeplitz structiore of the 

covariance matrix Tj^ allows (2.42) to be solved in a 
recursive manner as suggested by Levinson [ 47 ). 

Here it need be pointed out that the Autoregressive (AR) 
modelling approach and the prediction filtering approach lead 
to the same spectral estimator. It has been shown by [ 48 3 

that MEM spectral analysis is equivalent to fittng an AR model 
to the random pi'ccess. This brings out an interesting fact, 
that the respresentation of a random process by an AR model is 
that representation that exhibits the maximum entropy. The 
equivalence between the prediction filtering approach when 
the present value of a r8.ndom sample is estimated by a linear 
combination of some past values, and MEM has been shown by 
[ 393 * “^he MEM approach can also be viewed as an attempt to 
find a filter structure that whitens the observed time series 
[ 353 . The whitening filter transfer function H(f) in this 
case is given by, 

N~1 

H(f) == ^ “k ®xp{-j2itfk At) ( 2 . 43 ) 

The above results for IIHIM spectral analysis is also valid for 
a comp lex -valued time series [49). It was also assumed that 

I* 

in solving (2,42) the R(k)*B for k < H were known. Burg [363 
suggested a method of obtaining in (2.42) directly from 

the data samples without prior knowledge of the autoc ©variance 
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function. Based on this method, Anderson [ 50 ] and Haykin [ 49 ] 
have developed a fast recursive procedure for MEM spectral 
analysis of a real~valued and complex -valued time series 
respectively. A Brief summary of the recursive procedinre 
developed By Haykin [493 is given Below. 

Re cursive Brocedure for MEiM Spectral Analysis of a Complex 
Valued lime Series : 

Equation (2.42) corresponds to fitting , an AR model 
of order H*1 or finding a prediction error filter of length I. 
^N -1 output power of the prediction error filter. 

Rather than finding the N~1 coefficients of this error filter, 
the coefficients of a prediction error filter of length 2 are 
first found. The coefficients of the prediction error filter 
of length 3 are recursively evaluated from these and so on. 

The equations to Be solved then are, 

m 

I ^ ®mk ^j.o * j 0,1,. . . ,te 

and a . = 1 for any m 

(2.44) 

where m denotes the order of the AR model Being fitted and 

m +1 gives the length of the prediction error filter whose 

coefficients are Being found. Equation (2.44) is now solved 

recursively By, increasing the matrix dimension from m to m+1. 

The initial value for P when m - 0 is given By 

m 

^0 & «(“) = 1 X I I 


(2-45) 
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The corresponding prediction error filter is trivial with a 
single weight equal to unity. After solving (2.44) for 
(ni”l)th order of filter, the next step for m involves the 
determination of m+2 unknowns (E.(m), I’m) 

from the m+l equations (2*44). The additional condition 
required to solve for the unknowns is obtained by minimising 
the average output power or the circular error power with 
respect to This circular error power is obtained 

by averaging the power from forward filtering and that from 
backward filtering, i.e,, when the filter is conjugated, and 
the samples are run reversed in time. According to Haykin and 
Kesler [49 3 » this circular error power is given by 

■1 N~m m 2 

’■m = 2^ [ I ) I 

m 

+ |X(k+m) + I a X(k+m-j) I ^3. (2.46) 

j=l 

Minimising 71^^^ with respect to and writing down the 
results [49) we get 

H-m . 

-2 J, Bnilc °«k 

(2.47) 

J" 

where and are given by the following pair of 
recursive equations, 



3a 


A 

^mk ®(iii-l)(k+l) *^(111-1) (m-1) ® ( m-l ) ( k-fX ) 

The iteration starts with m = 1 and 


(2.48) 


®lk = 5C(k) 

011 ,= 


(2-49) 


The other coefficients of the (jn+l) long prediction error 
filter are calculated by 


%k ~ “(m-l)k “mm (m-l) 

(2.50) 

and the prediction error powers are given hy, 

^m= P(m-l) <2.51) 

The autocovariance function R(m) at each step m is obtained 

fey> 

R{m) *= « f a . R(m-k) (2.52) 

k=l 

The whole recursion is to be carried out upto and including a 
certain value of m « KAX. Proper choice of MAX. is an important 
part of MEM spectral analysis. Various criteria for the choice 
of MAX have been suggested [51-54 ). Ulrych and Bishop [42) 
observed that -spurious details were introduced into the spectrum 
when excessively large values of MAX were used. Small values of 
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MAX tend to smoothen out the resulting spectrum, thus 
annihilating some significant detail. However, the Jlnal 
Prediction Error ( EPE ) criterion introduced by Akaike [51“*53]I 
is the most widely used criterianto choose MAX, If are 

a set of AR coefficients based on a certain sample realisation 
of a time series, ilkaike’s EPS measures the error that is 
incurred on the average when these coefficients are used as 
the model for other sample realisations of the same time series. 
For zero~mean processes, the following estimate of has 

been developed by Akaike [51), 


FEE^ 

m 


N + ( m + 1 ) 

— p 

K - (m + 1) “ 


(2.55) 


The order of the AR model or MAX is chosen as one that minimises 
the FEE, A major objection to Akaike ’s EPS criterion while 
using Burg's recursive procedure is that the AR model orders 
tend to be overestimated when the range of search is close to 
the length of the time series. However, it has been observed 
to give good results [39,54) if the range of search is limited 
to N/2 where N is the length of the time series. 


Once MAX has been fixed, the autocovariance function 
R(k) for k > MAX can be extrapolated by 

MAX 


R(MAX + L) « - “MAX,k ^(MAX-k+L) ^ ^ 1 , 2 .,. ( 2 , 54 ) 

The inverse of the autocovariance matrix can also be found very 
easily as shown by Burg [jeJ* A^ be the NxN complex 
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a.utocovarianc e matrix of the time series of length 

"r( 0) R*(1) R*(2) ... R'^dT- 

R(l) R(0) R^(l) ... 


* • * « 

R(F-l) R(N>-2) R(T\T-3) ... F.(0) 

and 

R(p) = E[X(m) X**(m“p)] 

Let % te a NxR matrix whose columns are the prediction error 
coefficients obtained at each step of the recursion, 




( 2 . 56 ) 
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and let be an jJxH ' diagonal matrix whose elements are 
the prediction error powers 


N 


^M/a 


■M 






V 

0 


f 

« 


then it can be shown that 




or P-1 I*'' 


( 2 . 57 ) 


( 2 . 58 ) 

( 2 . 59 ) 


Regarding the statistical properties of the MEM estimator, 
due to the nonlinear nature of estimation, the variance of 
the spectral estimate cannot be easily obtained. Iromer,[553 
and Berk Ct63 determined the asymptotic properties of the 
MEM spectral estimator. According to their findings, MEM 
spectral estimate is Gaussian and unbiased in an asymptotic 
manner. If S(f) denotes the MEM estimate using a model of 
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order m and time series of length N, then, 

Var[§(f)3 = (2m/N) S®(f ) 
where S(f) is the true spectral density. 

Estimates for various moments of the MEIM spectral 
estimator are difficult to obtain because of its nonlinear 
nature . 
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OHiLPlCBR 5 

O PTIMUM HEGFJVER PORMULATION 


In this chapter y the optimum receiver for the radar 
detection problem posed in Chapter 2 , is obtained. We 
consider only Swerling 1 target models and the Feyman-Pearson 
criterion for optimisation. In Section 5.1, we derive the 
optimum receiver by solving an associated integral equation. 
In Section 5»2, we obtain an expression for the performance 
of such an optimum receiver. 


3.1 I ntegral Bquation Approach to the Optimum Receiver 


Restating for convenience the radar detection problem 
posed in (2.3<3), we have for a Swerling 1 target, 

M-1 

r(t) = yp. J S f(t~mT ) exp(jmjjt) 
m=o ^ 


M~1 L ^ 

+ I J \k ) + w(t) 

^ m=o k=o ^ ^ 



^ M~1 H, m, ~ 

r(t) = I ^ \lc f(t~lcT-mTp) + w(t) 

Ti < t < T^ (3.1) 

where B[f) « 0, E[^!r(t)3 * 0, * 0, B[|b 1^ ) = 2cr| , 

B[w(t)w’^u)3 « NQ6(t-u), ECbjjjj^b*^) iC(Di-n) &^y 
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Eq-uation (3.1) 
convenience, we 


ujses the notation of Section 2.4. 


denote (m«“n) as Rjj.(m-n), 
Ic 


For 


This hypothesis testing problem of detecting a target 
in the zero’ th range bin and with Doppler shift is a case 
of detecting a signal with unknown parameters in coloured 
noise. The usual way [ 5 3 to solve such problems is to use 
the Karhunen'-loeve (ICL) expansion r(t) = where the 

9 j^(t )*6 form a complete orthonormal set. 9 jj.(t)’s are obtained 
by solving the integral equation, 

\ ^ K~(t,u) (Pj^(u) du (3.2) 

^i 


where Z-(t,u) is the covariance function of the coloured 
noise. The optimum receiver computes 


T. 


A i / r(t) 


T. 

1 


g*(t) 


dt 


(3.3) 


and compares with a threshold. g(t) is obtained by 

solving the integral equation, 

M-1 ^f 

I f(t-mT ) exp( 3 Wj^t) = J 1fe(t,u) g(u) du ( 3 . 4 ) 
m^o ^ T^ ^ 

< t < T^ 

We now proceed to solve the integral equation (3*4) and 
obtain g(t). S^(t,u) is given by 
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K£*(t,u) 

where 




E[(H„{t) + W(t)(a*(u) + 5*(u))) 


M-1 I, . . 

I I t(t-kavml ) 

m= 0 k= o ^ 


(3-5) 


Ab thermal noise and clutter are independent processesj 


%(t,u) =s S[n^(t) n*(u)3 + E['C?(t) w*(u)) 


= K^(t,u) + NQ6(t-u) 


( 3 . 6 J 


where K (t,u) is given hy, 

w 


K^(t,u) = B[nQ(t)nJ(u )3 


M-1 I „ 

Kt-kO^-inTp)]^ 


M-1 L * * 

[ I I Ki ? (^~ 3 T-n.T )33 

nrO 3 «o ^ ' 


(5.7) 


As = R]^(m-n) 6^^^ from (3.1), we have, 


M-1 M-1 1 

e (t.u) = s+ y I I V“-“) f*(u-kT-nT ) 

m=o n?=o Is&o ir' 

(3.8) 

Prom (3*4) we have 




s(t) 


T / - 1 M -1 ^ 

“ ^ K (t,u) g(u) du + ^ I f(t- 


No "c 


^ I ii^r-mT )exp(ow t) 
^0 maO ^ 


T . < t < T ~ 

1 I 


(3*9) 
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This is a Fredholm integral equation of the second kind with 

given by ( 3 . 7 ) as its kernal. As K:^(t,u) is Hermitian 
and non- negative definite (being a covariance function), it 
can be shovm [57) that there exists a unique solution for 
g('t). 


Substituting (3«6) and (3.8) in (3*4) we get, 


M-1 

I ?(t-mT ) exp(jw^t) 

m=o ^ 

M-1 M-1 L , ^ 

= Ex M I I I %(in-n)f(t-kT-mT^)f*(u-kT-nT ))g(u)du 
^ Tj_ itt=o n=o k=o ^ E ^ ■ 

+ NQg(t) (3.10a) 

or 

M-1 ^ 

I f(t-mT ) exp(aw^t) 


M-1 1 m-1 ^f^* 

I I (f(t-kT-mT ))[ y Sj^(m-n) f t (u-kT-nT )g(u)du} 
m=soli£:o •^n=o ^ 

+ N^g(t) (3.10b) 


Observing the term within the square brackets in (3.10b) we 
find g(t) will be of the form 


^ ^ M-1 , 

g(t) = I f(t-mTp) exp(jWpt) 

Hts o 

M-1 I 

" io io 


( 3 . 11 ) 
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This gives the structure of the optimum receiver. ¥e will 
now substitute this expression for g(t) in the integra,! 
equation (^.lOh) and solve for and to determine the 

receiver completely. Define 


fj(t-kT) = 


and 


f(t - kT) 


f(t - kT - Tp) 


f(t - kT - (M-l)Tp) 


Mxl 


for k = 0, ...» L 


Sk(i) 


\( 0 ) 






Ejj<M-2) ... 

for k = 0» . 


Rk(o) 

, L 


(3.12a) 


Mxm 


(3.12b) 


Writing (3»8) iii matrix notation we have^ 


K 




Q(t,u) fJ(t-kT) Aj^ f£(u-kT) 


(5.13) 


Define, 
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f(t) exp(oui)t) 

' 

1 

fg(t) * 

■ 

f(t-Tp) exp(3wx)t) 

* 

• 

• 

I — 

, ig - 

1 

• 

« 


f(t-,(M-i)T ) exp(DWjjt) 


• 


and 


~ 4 


& 


ok 


^Ik 




Mxl 

(3.U) 


for k = 0, 1, 2 , , L (5.15) 


Mxl 


Writing (3.11) in matrix notation we have, 

g(t) = 643 + jJ, fjCt-M) 4 


(3*16) 


De fine 


A 


5j^g = J fg(t) dt 


1- 

■^1 


T. 


T- 

X 


f*(t-kT) 
f^Ct-kT-T ) 




# 

S* 


f*(t-k3>-(M-l)Tp) 


Jf(t)exp(jw^t) 


f(t-(M-i)a?p)x 


exp(jWpt) )dt 


for k s= 0,1, . . .1 


(3.17) 
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Now because of the conditions imposed on the •waveform ?(t) 
in ( 2. 29) will be a diagonal matrix. is given by 



where ( 3 * 18 ) 

^f 

^kS ~ i ^*('t'-M~wTp)f(t-wTp) exp(aa)jjt) dt 
•^i 



(3.20) 
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where is an identity inatrix and 

^k-h ~ ^ ?*(t~kT) f(t-hT)dt. 

^i 

¥e Observe = i and 

Substituting (3.6), (3.13) and (3.16) in the integral 
equation ( 3 . 4 ) we get. 


f^(t) I 


S 


Sic) 

T 


f L 


■•■hj [J. 


k£o 


^%s ■*“ I ^l(^~'li)ga!)du 

0 ^ 


or 


(3.21) 


♦s^ rp 

fs(t) I 


s 


N 


1 . 


s ' lio "I' 


^ io V "s 

+ ®t ^ J "^(t-M) A 5 g 
q =0 kso K Kq q 


( 3 . 22) 


he fine 


0 a P T 
^kS kS 


for k = 0,1, 2, . . . ,L 


(3.23) 



45 


From (’3.22) we have 

fg(t) Ig » Nq ?g(t) Is + I ?j(t-kT) g 

l<*o 
L 


+ ®t 5a, I £^(t-kT) c. 


k==o 
L L 


kS 


+ \ l l fi(t-kT) I -B 3 


4=0 k=0 


q "q 


(3.24) 


Observing (3.1l) we see that g is just a scale factor and 

as 

by choosing a particular value for g^g/ the and g(t) 

will be suitably scaled. In conformity with (3.10) we 
choose ‘ Substituting the valve of » l/N^^ 

in ( 3 . 24 ) we get 

L r EL 

° = ”o fl(t-.kT)-g^ + ^ i*(t-kT) C^g 


q=o k»o 

Expanding (3.25) we get 


* \ _L ,1 K 8k., 5, 


(3.25) 


f J(t) [StAo+'^o^D^So + ... E,. 5j^Sj(t.LT) 


+E^ 6l il(t) % 9i +il(t-T) [Ej. Aj+N„Ij9i+. • -St ^L®1 


-T 


T, 


ft* -"T / . « . 

+E^p3^fj(t) AnS 


.T 


o^L 


+ fj.(t.W)[E^Aj^+N^Ip]gj_ 




t r t-T 


-T 


N. 


[f (t) A Cq 5 + .. , f j(t»XiT) (3.26) 
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Now miiltiplying (3*25) or (3*26) by for i=0,l,2. ..L 

and integrating with respect to t from to we obtain the 
following set of matrix equations, 


I L 


f'o 6l-k % + 


Et I 


J \-k Ak ^ks for i=0,l,2,...L (3-27) 

0 Kr— 0 


Or rearranging termS] 


1 I 


E 


I 


t 5! I ^i~k ^k-q^k ®q ^i-k^^o ”*■ ®t %^®k 

qs :0 k=o ^ ^ k~o 




B X) 

“ “ J ^i-k \ \s for i=0,l, 2.. .h (3»28) 
0 1&:0 


Again rearranging terms and expanding, 


Is/o 


■'' l^Pi-l^VD **■ ^^t ^1^ ®t^i-k^k-X "^k^^L 

Ks= O 


k^l 


I^^i-X^Vp + ®t '*■ Jo ®t^i-k^k-*L 

kj^L 


®t ^ - 

^ Jo 


for i e 0,1, 2* . .X 


(3.29) 
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Equations (3.28) or (3.29) are a set of (L+i) simultaneous 
matrix equations in (L+i) unknowns for k = 0,1,2,...L. 

The g^'s are obtained as solution to this set of equations. 
We observe that C^g is the only term that will change if 
the hypothesis testing problem is changed to detect targets 
at other range bins. 


In an actual situation, the number of matrix ■'iequat ions 
to be solved depends on the €^^ 3 , and terms for k=0, 1 . . ..L. 
From (3.18), (3.2o) and (3.23) we observe that 3^^ and C ^5 

terms depends on the ambiguity function 


I / f(t 

T. 

1 

of the waveform f (t) 
shown in Fig. 3.1. 


- ‘5^'^*(t + Jl)exp(j w't)dtp 
2 2 

. The optimum receiver structure 


is 


3 . 2 Performance of the Optimum Receiver 

From Section 3.1/ the optimum receiver computes 

Aw/ r(t) g*(t) dt and then compares UP with a 
T. ' ' 

s , . X . 

threshold, 




f(t) g*(t) 


dt 



r\ 


T 


(3.30) 
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* 1^-1^ 

§ (t) = -Y“ I f (t-mT )exp(-3a3j.t) 

0 ms=o ^ ^ 

L ^ ^ 

+ Z Z f (t-kT-mT-) 

ra^o l&=o P 


Fig. 5*1 Optimum Receiver for Swerling 1 

Target Models (Complex Operations). 


ffiA V 
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The sufficient statistic ^ is a complex SauBsian random 
■variable. jl|^-will have an exponential probability density 
function. In such a case, it can be shown [28 3 that 
1+ A 

= (Sd) (3-31) 

where Pp denotes the probability of false alarm, and 
Pp denotes the probability of detection, and 


E[ U| 2 /H ) - E[ \l\^ /H^3 


/Ho5 

Now 


(3.52) 


B[| /Hq) = Et / f r(t)r*(u) ~g*(t) g(u) dt du] 

Ti Ti 




T 

J 

T.- 


f 

i^(t,u) 


g*(t) g(u) dt du 


(3*53) 


Using ( 5 * 4 ) in (3«33) w© gei» 


"^f M~1 . 

1^ 1^ = I [ I ) exp(j wj.t)3g*(t)dt 

T^ m= 0 . ^ 

% M-1 * 

= I [ I i'(t-niT )exp(jWj.t)3 g(t)dt (3.34) 

Tj_ ms=o ^ ^ 


Also, 



M -1 


j 0 


M -1 


^ ^ (^-5aT )exp(-ja) u)]g(u)du dt 

xte 0 ^ 


B' 


+ E[|J |2 /H^} 


(3.35) 


Prom (3.32), ( 3 . 34 ) and (3.35) «e get 
(P 

2 f M ~1 

A« 2ag f(t~inTp) exp(ju)jjt)3g*( t) dt 


(5-36) 


substituting for i(t) from (3.II) (with have, 


T. 


f M -1 

' ^‘'b ®t J C J, f(t-mT^) expCjMjjt) )x 


111=0 


1 M-1 

”E 


M-1 Jj 


0 Jo ^ l ^ f *( t - kao - pTp )) 


Simplifying, 

A 3 . orT® p r M , M“1 ^ m , 

20 b E, 

1 


(3*37) 


(3*38) 


m 7f - 

where f *( t-kT-mTp)f( t-mT^) expaw^jt) dt 
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The degradation in performance when clutter is present is 


due to the pjjg terms. Of course, also depends on 
tihe terms as can be seen from (3.29) . As is evident from 
( 3 . 29 ) and ( 3 , 39 ), the optimum receiver performance depends 
on the statistics of clutter and thermal noise as well as 
on the ambiguity function of the waveform f(t). This 
performance of the optimum receiver could be improved by 
waveform or signal design. Of course, this calls for 
1oint optimisation of the transmitter and receiver for a 
particular clutter environment. 


As in general, the statistics of clutter are time 
varying, the coefficients (g^’s) ‘of the optimum receiver 
have- to be continually updated to meet this changing clutter 
environment. This makes implementation of the optimum 
receiver unwieldy and complicated* If there are 16 hits /sc an 
(M=16) , a set of (L+i) 16 simultaneous equations with complex 
coefficients have to be solved each time the receiver 
coefficients (g^*s) are to be updated. Also in practical 
situations, where the statistics of clutter are a priori 
not known, optimum receivers which continually 
estimate the required clutter statistics and then 
update their coefficients (g^'s) are difficult to 
implement with present day technology. In view of this, 
in Chapters 4 and 5, we propose and evaluate suboptimum 
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receivers and study how the J^axxmum Entropy m.ethod (ffiM) 

can he used to adaptively determine the parameters of these 
receivers. 
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OHAPGSlfi 4 

S TRUCTURED OPTIMUM BECEIW.R gORMULATlOII 


In Chapter 3> "WQ fonnd that the optimism receiver 
implementation in a changing clutter environment is 
extremely complicated. Global adaptive optimum receivers 
which continually estimate the required clutter statistics 
are difficult to implement with present day technology. 

In this chapter, we impose a structure on the receiver and 
within the constraints imposed by that structure, derive 
an optimum receiver for Swerling 1 target models. This 
receiver is termed Structured Optimum. Receiver. In 
addition, we see how the scalar complex form of the 
maximimi entropy method (tIElM) can be successfully 
applied to continually estimate the required statistics of 
clutter therby enabling the parameters of the structured 
optimum receiver to bo determined adaptively. In Section 4.1, 
we discuss the structure that is imposed on the receiver. 

In Section 4.2 we obtain a structured optimum receiver. for 
Swerling 1 target models and compare the method followed to 
those of Delong and Hofstetter [Id], Spafford [173 and. 

Rmmler [163* In Section 4.3» we obtain an expression for 
the performance of such a receiver. In Section 4.4, we 
evaluate the performance of such a receiver for the cases 
when the required clutter statistics are a priori known and 
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when MSM is used to estimate the required clutter 
statistics. The performance of the structured optimum 
receiver in hoth these cases is compared with that of a 
conventional matched filter receiver. A large numher of 
receiver operating characteristic (ROC) graphs indicating 
the performance of such a receiver, under varying conditions 
of the target, clutter and thermal noise are included in 
this section. These ROC's have been obtained by computer 
simulation and numerical computation. 

4*1 Matched Filter Structure Imposition 

Restating for convenience the radar detection problem 
posed in ( 2 . 30 ), we have 

M-1 

r(t) = I f(t-mT ) exp(jwpt) 

ntso ^ 


M-1 L 



n(t) w I I \ir :f(t~kT-mT ) + w(t) tK 

m?=o k»:o ^ ^ 

T^< t (4.1) 

where E[b 3 = 0, E[w(t)) = 0, « 0, E[|^|^3 = , 

E[w(t) wtu)3 = Eq 6(t-u) , = R^ (m-n) . 

Jx 

Equation (4*1) uses the notation of Section 2.4. Eor 



55 


convenience, we denote (m-n) ae R-(m-n). 

°k ^ 

The received signal r(t) is passed through a matched 
filter (matched to a single pulse return) and sampled 
appropriately for returns corresponding to a single range 
bin (aero’th range bin in this case) during the entire scan. 
This is shown in Pig. 4.1. A matched filter structure is • 
chosen, as intuitively it appears to be the best structure 
to impose and relatively easy to implement. Also a matched 
filter is the optimum receiver in the absence of clutter. 

It turns out that in the case of a Swerling 1 target model, 
imposing this structure implies choosing a class of transmit/ 
r e ce i ve wave f orms . 


r(t) 


MATCHED 
FILTER ■ 

Mr 


OPTIMUM 
.PROOF. SSOR 

■ — 

OH SAMPLES 

Mt) 

m 

^ A 


h(t) « f *(TQ-t)x Sample at 

®xp(~jWjj(T^-.t) ) tsT^ + mTp instants} 

m = 0,1, . . , ,M--1 




5l 


n 


T 


Pig. 4*1 Matched Filter Structured Optimum 
Receiver (Complex Operations). 
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IhG impiilse response |i(t) of the matched filter is given by 
^(t) * f (T^-t) expfj6jjj(T^-t)) (4,2) 


The output Of the matched filter r at time t 


m 


on hypothesis ie given "by 

f M~1 

^m == Bjj^flu-mTp) exp(3a)j)U)3 




p 


[ f*'( )3'i^ 




■f M-1 1 ^ 

+ j/ t'®t J, J, \k 


t«:T ^-mT 

Jr 




mwo kso 


[f (Vt+u) exp(-3w^(!rQ-t+u)),'|du 


+ ^| w(u)i*(To‘-t+u) exp(-ja>j5(To-t+u))du 


t=!r^4-mT 

O p 


t=T +mT 
0 ‘ 


(4.5) 


Using the conditions stated in (2.29) and simplifying 
wo haVG » ■■: ; ■ ■ , , . . ,.. 


■m 


^'t ®m ®2cp(:i«^pni®p) 




jJo ^mk J V--niTjj)exp(-.;3oj^(u-ml ))du 


t W 


m 


or 


P' 
(4.4) 
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T. 


^mk I f(u-kT)f*(u) exp(-ja)pU) du 


T. 

a 


+ w. 


where 


m 


T. 


%* ^(■t) ^ (t-mT^) exp(-jo^(t~mT ))dt 

i 


(4.5) 


(4.6) 


De fine 


~ ?(t-.kT)r(t) exp(-jwpt) dt 


T. 


and 


d^ ss 
Ql 


Jo 


(4.7) 


(4.8) 


It iB observed that depends on the ambiguity function 


Of the waveform f(t). 


Prom (4.5) » (4.7) and (4.8) the hypothesis testing 
problem can now be stated as 


®’'P(3“e'“®o) + V^t aii t : E 


'm 


VEt * «. 


m * ^0 

m * 0,1,... ,M~l 


(4.9) 
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where b[bJ * 0, B[d„) » 0, = 0, E[|Sj = 2o-=, 



0 mn 


The above hypothesis testing problem is now solved to 
obtain the structured optimum receiver* The Neyman-Pearson 
criterion is used for optimum processing of the samples. The 
receiver structure so obtained will be a structured optimum 
one as r(t) is first passed through a matched filter before 
solving the hypothesis testing problem. 


4*2 Structured Ontlmum Receiver for Swerling 1 Target Models 
Prom (4*9) for a Swerling 1 target model we have, 

” ^t '^m * 

m = 0,1,2,.. . ,M-1 (4.10) 

The are jointly complex Gaussian random variables as 

they have been obtained from a complex Gaussian random 
process, let 

Using (4.11) in ( 4 . 10) we have, 

* ^t ® ®xp(jw^mTp) + : % 

w. ni 0 


(4.12) 
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Ijq’s are also jointly complex Gaussian 
Define i; 


random variables. 


^f. 


f 4 


■M-1 


f d ^ 



• 




^0 




• 

> W 

• 

• 


• 

• 

a. 


mt 

L 





«<« 




. 

A 

a i 

U 

^1 

* 

♦ 


1 ' 

• 

• 

' ■ 

- A 
f r = 

^0 

IK# 

• 

« 

How 

_V1 


®xP(Jwjj(M-l)3^) 


‘ 

/m~1 


B[d3 = 0, B[i3 = 0, E[;3 = 0, E[B] = 0, E 

Let 

®[4 d**3 = ah* I 


( 4 . 13 ) 


E[w3 * 0 (4.14) 


(4.15) 


where the term 2h= is used to make the diagonal elements of 
■■ qual to unity. 2h represents the expected value of the 
clutter power in each sample when E . = 1, 

X 

Ifow 


and 


I = 1®+ d + w 


B[I I ^3 _ B[d 1*^3 +E[w w*^3 E . 2h* I 


(4.16) 


where ie an identity matrix. 


■d " "0 
(4.17) 
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l©t 

E[f 1*^3 » Aj (4.18) 

where 2P ^ * 2h“ and tise diagoml elements of Ij 

are equal to unity. 2P * represents the expected value of the 
total clutter and thermal noise power in each sample, Aj is 
a positive definite matrix. 

Using (4.12), (4.15), (4.15)» (4.16) and (4.18), the 
hypothesis testing problem can now be stated as 

r a B u + I ; 

(4.19) 

r * I : 

where the probability density function of f is given by 


p-(f) 1 exp[- -i- i)3 (4.20) 

(2xp»)“Det(Aj) ^ . 


Ihis hypothesis testing problsn will be solved using the 
Neyman-Pearson criterion for optimisation. As stated in 
Chapter 2, this criterion leads to the likelihood ratio test. 


A(?) 




< 

H. 


T), 


(4.21) 


where /H joint probability 

density functions of the observations given that and are 
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true respectively, A (r) is the likelihood ratio. We first 
obtain the likelihood ratio given B (A(?/B))f Now, 


A(f/B) 


where 




(4. 22) 




Bet(Aj) 


exp[*’ ■i— -(r- ySx B u) ^ X^(r“V^j, S u)3 

2^2 V A T 

(4.23) 


and 




- ej5,[- ’’ 5)7 

A \ r>%n 2 A 


2 P 


(4.24) 


Substituting the values of and Pp/jj (?/Hq) in 

(4.22) we get, 

A(r7B) = exp[— bV'' r 


2 P^ 

^ ni -—1 ^ 

•n » .i- A -^s 


+ B A”-^ u - B^lBj * u*^ Aj^. u)) (4,25) 

As is Hermitian, the quadratic form u*® Aj^ u will be a 
real number. let 

T « u*^ u 


As ? ^ u and u*^ r are scalars, 


(4,26) 
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u)* (4.27) 

u) + Sf 

|»|‘I)3 (4.28) 

Let 

B « I B| Gxp(je ) (4.29) 

#» ^ 

As B Is a zero mean complex {Jauesian random variable ^ )b1 

has a Eayleigh probability density given "by 

; I S 1^ 0 

(4.50) 

; |B! 5.0 

0 is viniformly dis-tributed between o 2it and is independent 
of |b|. let 

u I exp(-ja) (4.31) 

Then 

a( 5’/I) =A (V|8| . 8) 

= e 2 p[ I 9xp(j e) I a- u| exi>(-ja) 

2p ^ ■ 

+ P^\i\ exp(- 3 e) j ulexp( 3 a)~B^jB j^Y)] 

(4.32) 



^*T 7*“1 ~ I x***l ^■vT / *"*1 
u Aj. r *= (u * Aj r) *s (r 

From (4. 25) 

A (f/S) = bxs[-A-(P. Btp*" 
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or 


A(r/B) = Qxp[ u| cos(0*a) 


- B. |Bl* y)) 


Now 


A(Ty|B| ) 


2tc 


2 % 


i A (2^/!B j ,9) de 


(4.33) 


(4.34) 


or 


- T 2it 1 B I* E. Y 

A(r/|B! ) « J exp[-- - 7 ^ 3exp[ 


2^ 


Aj^u| cos(e -a) 




2P 


P 


} 


Where is the modified Bessel function of zero order. 

Now 


A(^) « ^ ir/\B\) Pj^j(|B|)d|B| 


(4.36) 


Substituting for P|;g|( lB|) from ( 4 . 3 O) we get; 

^ 1 = 1% . Ir 1^1 

“J io( : ) X 


A(r) a= J exp[- 
0 


2P 


3 


exp(- dlBl 


2d 


( 4 . 37 ) 


B 


This is a special case of the integral, 

« - exp[~3,B2( ^ +1; v+1; •—) ( 4 . 

2 p''r(v+l) 4 p®' ^ 4p^ 


38) 
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where r(.) is the gamma function and e; a) is the 

confluent hype rgeome trio function. 


Simplifying (4,57) using (4,56), and taking the 

log likelihood ratio, we obtain the optimum decision riole 
as 




let 






^ "xT . Ti-'l •* 


(4,40) 


(4*41) 


The optimum decision rule 


T ,7 .a 

% 


r V p 

> 

< 

«0 



can then be written 


or 




h IP 


h<3.uation (4.45) can also be written as 


(4.42) 


(4.45) 


M-1 % 

I f ^ ^'^12 ^ 

'uteo ^ (4*44) 

o 

The structured optimum receiver corresponding to (4,44) 
is shown in Pig. 4.2* 

We observe from (4.5) that could also be obtained 
by correlating r(t) with f*(t-inTp) exp(-«j(jjj,( t-mT )), 



Matched . 1 

^Filter. : 


Sample at 
t =s 4- ml 


TiiPlBD DELAY III® 


Close at 
t«!Po+(M«l)T 



"1 
> - 
< 

Ho' ■ 


V* C T1, 
mm' p. , 


Eig«4. 2 Str-uctiired Optimum Receiver for Swerllng 1 
Target Models (Complex ..Operations). 


V 
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(4.45) 


Using ( 4 . 45 ) in ( 4 . 44 ), the optimum test is given by 


„ M-l 

I / r(t)[ I V exp(3ai mT ) f*(t-niT ) x 

HfcrO . -P -P 


2 

®xp(--jWjjt)3dt I < Tij (4.46) 

An alternate realisation of the receiver can be obtained by- 
correlating the received waveform r(t) with 


M-1 

I V expCowr^mT ) f*(t-m!P ) exp(-3ajT)t) 

10=0 " ■t' ^ 

and taking the modulus square to get the sufficient statistic. 
Pig, 4.3 shows an alternate realisation of the structured 
optimum receiver based on (4.46). 

From ( 4 . 46 ) we observe that imposing a matched filter 
structure on the receiver is eqtiivalent to choosing a class 
of transmit /receive waveforms. In this particular case, the 
transmit waveform (by transmit waveform, we mean the reflected 
signal from the target only) is given by 

M-1 

exp(j(opmTp) f(t-mTp) exp(ju)j5( t-mf^) ) 
and the receive waveform is gi-ven by 


(4.47) 




Vq f (t) exp(3u>jjt) 


‘ exp(j o^(t~Tp)) 



d-t 

T. 

1 


Vm.i f*(t-(M-l)Tj,) X 

exp(d<^(t-(M-l)Tp)) 


Fig. 4. 3 Altermije Realisation of the Structured Optimum 
Receiver (Complex Operations). 
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M-1 

'I f(t-mTp) exp(Ja^{t~mT )) (4.48) 

IXt?* 0 

Also it need be mentioned here that Delong and Hofstettor [18^, 
Spafford Cl7}> and Rummler [16) have obtained the same receiver 
structure for a slowly fluctuating point target under a 
different criterion. They have assumed a class of transmit/ 
receive waveforms similar to the ones in (4,47) ^nd (4*48) and 
maximised the signal to interference ratio (SIR). Bach has 
followed a different method of maximising the SIR, which is 
their criterion, and the results obtained by them are identical 
to the one obtained here. It can be easily shown that when 
the return signal r(t) is a complex Gaussian process, 
maximising the SIR is equivalent to maximising the probability 
of detection (P^) for a given probability of false alarm (Pji). 
Delong ot*al[ island Spafford [173 have also proposed steps 
for joint optimisation of the transmitted signal and receiver 
for a class of transmit/receive waveform for a given clutter 
scattering function. 

A knowledge of v or Aj is required to implement 

the receiver shown in Fig. 4*2. The inverse covariance matrix 

~-l ~ 

Aj may be estimated from using the complex form of the 

maximum entropy method, (MEM), during tho<^e scans, when no 

target is present in that range bin, or H decision to that 

effect has been taken. It was shown in Chapter 2, that 
’“*•1 

A]- could be obtained directly without the need to obtain 
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*r « 

and then compute its inverse. This estimate of could he 
averaged over several such scans. In a recent paper [45)» 
Koslur and Hay-kin report using recorded radar clutter data 
to carry out MEM analysis to obtain the clutter spectrum. 

It was found that 16 samples (or equivalently 16 hits/scan) 
were adequate to get a ’very good' spectrum. The difference 
in the spectrum, when 256 samples and 16 samples were used 
respectively, was negligible. Even 8 samples gave a 

mit T 

reasonably good spectrum. Thus a good estimate of may 
be obtained from real radar data by using MEM. A$ is 

being calculated during those scans when no target is present 
or a decision to that effect has been taken, the quasi time 
invariant nature of the statistics of clutter are also taken 
into account. Thiis the structured optimum receiver for 
Swerling 1 target models can be implemented even when the 
clutter statistics are a priori not known, 

4.3 Performance of Structured Optimum Receiver for 
Swerling 1 Target Models 

Prom (4.44) the optimum decision rule is given by 

i M-l ?1 

=1 I 1“ i It 

IQssO Hq 

where J is a sufficient statistic. 

5 Now are j ointly complex Gaussian random variables 

as Ijn’® jointly complex Gaussian and Ijj^’ s and B are 
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independent. il ls a linear combination of and is also 

a complex Gaussian random variable. 


The probability of false alarm Pj, --is given by 

CO 

Bp = f (2/Ho)d2 (4.50) 

llj 0 

and the probability of deteotion Pjj is given by 

oo 

Pp *1 Pg/H (4»5i) 

riT 1 


The probability density function of z vill be exponential. 
Then from (3.51) and (3«52) we have. 


1+ A 

Pp = (Pp) 


(4.52) 


where 


E[ I'M^ - E[ |t|‘' 
E[ br /Ho] 


(4.53) 


A is called the signal to interference ratio (SIR)* 
It may be observed from (4.52) that maximising P^j for a 
particular value of P^, is equivalent to maximising A* 


In the case of the structured optimum receiver. 


M-1 . 

2 - I 


ltt=0 


m m 


HfcsO 10*0 


**T 


(4.54) 

(4.55) 


V 
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E[ 1“ /Hj,] = J r = V*® B[i f ®]v 

= V Aj V (4.56) 

E[|4|" /%) = B[(v*® r I*® v)/H^O = ^?® E[{r f ®)/H3^3v 

(4.57) 

Now 

E[r5* = E^ E[ |B p u u*^ + 2p \ (4.58) 

Jj& t; 

A - ~ 

Ag, * u u 

Then 

[E[ lil V%) = 2 a| ?-® Aj f + v'*® Aj v 

and 

2c7| V 20^ V*® Aj v 

~ 2 P» v‘® 2i V " 

2o’g B^ and 2P ® are the expected value of the target and 
interference power in each sample respectively. 

4*4 Performance Evaluation of the Structured Optimum Receiver 
by Simulation and Numerical Oomputation 

The performance of the structured optimum receiver, 
for both the cases when the reciuired clutter statistics are a 
priori known, and when the maximum entropy method is used 
to estimate the required clutter statistics, is evaluated by 


(4.59) 

(4.60) 

(4.61) 
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means of computer simulation and numerical computation# The 
performance of the structured optimum receiver in both these 
cases is compared with the performance of a conventional 
matched filter receiver# Performance evaluation is done for 
varying conditions of target power ^target velocity# clutter 
power# clutter spectral width# thermal noise power etc. As 
is commonly done in most performance evaluations# the clutter 

returns are assumed to have a Gaussian shaped power spectral | 

I 

density# A Gaussian shaped power spectral density corresponds ! 

! 

to an infinite step autoregressive model. j 

The entire simulation scheme is discussed in Appendix A. 
Briefly, zero-mean complex Gaussian clutter samples dj^'s, with 
Gaussian shaped power spectral density, are generated on the | 
computer. Zero-mean complex white Gaussian samples Wj^.’ s are, 

1 : 

also generated on the computer and these are added to djj.'e to [ 

•r •• i: 

obtain interference samples These Ij^’s correspond to I 

the output of the matched filter in Fig. 4« 2 when no target | 

is present. Then the complex form of the maximum entropy | 

method is used to estimate the inverse of the total interference | 

t 

fv 

covariance matrix. Performance • evaluation for various cases | 

is now done using variations of (4.61). Details of this ! 

are given below. 

Number of hits/scan (M) *=16 

Pulse repetition frequency (l/Tp) ts 54 O Hz 
This means each sample is spaced 1/540 seconds apart. 
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We use tiie following notation: 
denotes the actual interference covariance matrix. 

'“—1 

denotes the MSM estimate of the inverse of the 
interference covariance matrix. 

denotes the expected value of the clutter p over /sample . 



Pj, denotes the expected value of the target pow3r/sample. 

^IT “ ®t* 

Nq denotes the expected value of the thermal noise 
power /sample . 

f^ denotes the two-sided half-power width (in Hz.) of the 
clutter power spectral density. 

f^ denotes the centre frequency (in Hz) of the clutter 
power spectral density. 

fjj denotes the Doppler frequency shift (in H?) due to the 
target motion. 

Aj and u are as defined in ( 4 * 13 ) aiid (4.59) respectivsly. 


The performance of the structured optimum receiver, 


when the required clutter statistics are a priori known, is 
evaluated hy calculating where 


m ^ ^ 

20 ^ E. v*-*- A™ V 


3 -"t 


m ^ ^ 

I’j A^ V 


(2h'B^.+N„)v*%v (Pe+Ho)^'"' Jj T 


(4.62) 
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wiiGre 


■*-1 *“ 


P 




1+A 


th 


(4.65) 


Por any given clutter statistics, this performance is the 
maximum attainable by the structured optimum receiver. This 
receiver is denoted as structured optimum receiver SI. 

The performance of the structured optimum receiver, 
when the required clutter statistics are estimated using the 
complex form of the HEM is evaluated by calculating 
where 


A 


2cr^ E^. V V 


" '(2Ji“E^+K^)V® Ij. t 
whore t = u 


^ m A# 

V* 4^ V 
(Pc+Nq) v*'^ Ij. V 


(4.64) 


Ppc: 


(I’d) 


1+ A 


stb 


(4.65) 


This receiver is denoted as structured optim-um receiver 02. 

The performance of the conventional matched filter receiver 
is evaluated by calculating 

’*0} 

’ Am w. rm a Am U 

( 4 . 66 ) 


"mf 


2d| Am U P^ Am U 


"T 


T 




(2h^E^+NQ)u*^ Aj u 


(Pc + Nc)ti^ Aj tt 


P 


(?])) 


1 + 


^f 


(4.67) 


performance of the conventional matched filter receiver 
no clutter is present is evaluated by calculating 
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^ideal * 


9 "♦T’ " 

T •• 

(4.68) 

Ho u 

No u 



(4.69) 

‘sto’ Amf 

■^^deal calculated 

for a largie 

cases by varying Pj, P^, N^, f^^ and f^. 

These are 


presented in Table 1. The receiver operating characteristic 
(ROO) graphs for some of these cases are given in Figs. 4. 4 to 
4.16. ' 


The following observations .are made on the overall 
performance of the structured optimum receiveri for both 
the cases when the clutter statistics are apriori known, 
and when they are estimated using the MBMs 

(a) In the fairly wide ran^ of situations considered 
the performance of the structured optimum receiver 
when cluttor statistics are a priori known or when 
they a 3 :^ estimated by MEIM, is far superior to that 
of the conventional matched filter receiver whose 
design ignores the presence of clutter in the 
received signal. 
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(b) In a large number of cases» the difference between 
^Bto ^th under 15 percent. This percentage 
difference decreases with increase in the target 
■velocity. MEM appears to be able to estimate the 
required clutter statistics ’reasona.bly well' -with 
only 16 samples. 

(c) The performance of the three receivers is greatly 
dependent on the target velocity or Doppler shift. 

^ideal 1 gives the maximum attainable 

performance when no clutter is present, everything 
else remaining the same. At very low target doppler 
shifts (20 Hz and below), the performance of all the 
three receivers is poor. However, even here the 
structured optimiom receivers perform better than the 
conventional matched filter receiver. At target 
Doppler shifts of 100 Hz or more and for f^ equal to 
10 Hz or less, the performance of the structured 
optimum receivers approach the maxim'um attainable 
performance (when no clutter is present) to within 
about 20 percent (in terms of A) even when the clutter 
is as high as 50 dB above the thermal noise power level. 
However, the conventional matched filter performance 

is very poor. 



The performance of the structured optimum receivers 
improve as the clutter spectral width decreases. This 
improvement in performance (in terms of A ) is very 
noticeable at target Doppler shifts of 40 Hz and below. 
At a target Doppler shift Of 40 Hz and f^ equal to 1 Hzi 
the difference in performance of the structured optimum 
receivers and the maximum attainable performance (when 
no clutter is present) is less than 20 percent even 
when the clutter is as high as 30 dB above the thermal 
noise power ♦ 

At a target doppler shift of 100 Hz and f^ equal to 
10 Hz, the variation in performance (in terms of A) 
of the structured optimum receivers, as the clutter 
power increases from 0 dB to 50 dB above the thermal 
noise power, is less than 20 percent. This variation 
in performance decreases as the clutter spectral 
width decreases. 



Pulse repetition frequency = 540 Hz ~ 

Number of hits/scan = 16 where a = Signal to interference ratio 

Centre frequency of the clutter spectra (f„) = 0 Hz 
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CHAPTER 5 



In this cha,pt^r, we propose a structured recoiver, 
that in raany instances is simpler to implement than the 
structured optimum rocoivcr of Chapter 4. This structured 
receiver is biased on interpreting the mo-ximum entropy 
nothod as an attempt to do terrains the structure of a 
whitoning filter for a time series* Again a matched 
filter structure is iraposed on the receiver for the 
reasons stated in Chapter 4. The received signal 
is passed through this matched filter (matched to a 
single pulse return) and the output of the matched filter 
is sampled appropriately for returns corresponding to the 
range bin of interest during a particular sco.n, Further 
processing on the resulting samples is done by a whitening 
filter. This whitening filter is adaptively estimated 
by the maximum entropy method from these samples at the 
output of the matched filter during scans when' no target 
is present or a decision to that effect has been taken. 

Of course, the performance of this receiver will be inferior 
or at most equal to that of the structured optimum receiver. 
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In Section 5»1» we prc^ose a structured receiver 
based on the above whitening filter concept and disctise how 
this structured receiver cen be used for both coherent and 
noncoherent detection. In Section 5»2 we obtain an 
expression for the performance of the proposed structured 
receiver with coherent detect ion ..for Swerling 1 target 
models. In Section 5.5> we present a generalised procedure 
and derive formulae to obtain the probability of detection 
(Pjj) and probability of false alarm (Pj,) for receivers 
where correlated complex Gaussian samples are modulus 
squared and added. This procedure also takes into account 
the effect of a whitening filter or a moving target 
indicator (MTI.) filter prior to modulus squaring operation. 
Using this procedure, the and P^ of the preposed structured 
receiver with noncoherent detection may be calculated. In 
Section 5*4i we evaluate the performance of such receivers, 
for both coherent and noncoherent detection, for Swerling 1 
target models. Their performance is compared with that of a 
double line canceller MTI receiver. This is considered 
appropriate as the structured receivers using noncoherent 
detection turn- out to be similar to an adaptive MTI receiver. 
Again a large number of receiver operating characteristic 
(ROC) graphs indicating the performance of such a receiver, 
under varying conditions of the target, clutter and thernal 
noise are included in this Section, These ROC's have been 
obtained by computer simulation and numerical computation. 
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5.1 Structured Reoei-ygr Based on the Whitening filter 
Oonceut of MEIM 

The received tiignal f(t) in (2.50) (4*1) is passed 

through a matched filter (matched to a sin^e pulse return) 
and the output of the matched filter is sanpled appropriately 
for returns corresponding to a single range bin (zero’ th 
range bin in this case) during a particular scan, !From (4*12) 
the resulting hypothesis testing problem for a Swerling 1 
target model can be written as 

Tjjj ySlt B exp(jWjjmTp) + 

^ m m o 

m » 0,1,. . . ,M-1 (5.1) 

The Ijjj'b constitute a coloured noise time series. As seen 
in Chapter 2, the maximum entropy method can be viewed as 
an attempt to determine a filter structure that whitens 
the observed time series. It can also be interpreted as 
an attempt to determine a prediction error filter for a 
given time series. Prom (2. 45)» "^^e transfer function H(f) 
of such a whitening or prediction error filter is given by 

(in complex envelope notation) 

H(f) * I exp(-D27tfiTp) (5.2) 

is:0 

Where a == 1. a.'s f or i *= 1 , 2 ,...>,MAX and MAX are determined 
hy solving )(2.44)rocursiVoly. In the hypothesis testing 
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problem of KAX are determined from rjjj'e 

during those scans when there is no target Or a decision 
to that effet ha,s been taken. This is illustrated in 
Pig, 5*1 where is the whitened output of the filter when 


m/^o 

WHITENING PILTEIR 

MAX 


H(f)= 1+ 1 a. exp(-j2TtfiT ) 

ittl ^ ^ 






Pig. 5*1 Interference Whitening Pilter 

Obtained by MEM (Complex Operations). 

the input consists of clutter and thermal noise. ITow 
whether the are actually white noise samples or not 

•ft* 

depends on the statistics of Ijjj*s and how adeqmtely an 
autoregressive (AR) model of a finite order can be used to 
describe the observed time series. In the rest of the 
discussion throughout this eectlon, e^^'s are assumed to 
be white noise samples in the context of the MEM whitening 
filter. The normalised covariance matrix of is 

therefore assumed to be an identity matrix. However, in 
Sections 5.2 - 5*4f when we evaluate the performance of 
structured receivers, the whitening filter is taken to be 
just a linear filter following the matched filter. In the 
current scan, when hypxthesis or is being tested, the 
r 's are sent through this whitening filter. Let x denote 

JJJ iU 
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the output of the whitening filter. During the initialisation, 
phase, to obtain for'm= 0, 1, (MAX-l), the past 

values of b for m = -1, -2, ... -{MAX), axe requiied. As 
these are not known, two variations of the hypothesis 
testing problem follow. The first variation is to assume 
the past values of r^j^’e form = -1,-2, .. .-(MAX) to be equal 
to their mean or expected value This is equa,l to 

zero as the samples are zero mean random variables. Then the 
hypothesis testing problem becomes, 

MAX - - 

*10 “ “i B ®*P(3 “i)(“>-l)Tp) + Ej, 

i^cin 

*m = '=m 

where cc^ » 1 and m = 0,1,..,, M-l 

This hypothesis testing problem will be denoted as hypothesis 
testing problem A, 

The second variation is to ignore the first (MAX-1) 
samples at the output of the whitening filter and test for 
hypothesis % or only from m = MAX, MAX+1,..., M-r-l, This 
calls for no assumption, but the entire available data 
containing target information is not being utilised. This 
hypothesis testing problem can be stated as 






(5.3) 
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MAX ^ ^ 

\ ^ V^'t ® ®xp(d Wjj(m-i)Tp) + :Hj_ 

X„ = c * FT 

® m * “o 

where = 1 and m =r MAX, MAX+1, . . . ,M-1 (5.4) 

This hypothesis testing problem will be denoted as hypothesis 
testing problem B, The difference in porformanco of atructurod 
receivers based on hypotheses testing problems A and B 
re spectively , is brought out clearly in Section 5.4. 


The hypotheses testing problems of (5.3) and (5.4) 
can now be easily solved. Ve will consider the structured 
receiver for the following cases: 


( ^ ) Evpothesls testinfc' problem A with coherent detection . 
Hypothesis testing problem A is solved using the Neyman- 
Pearson criterion for optimisation. A similar hypothesis 
testing problem has been solved in Section 4.2. Using the 
same procedure the decision rule can be stated as 


M-1 

1 I ^ I 


m =:0 
M, 


< 


'T 


where s, 


m 


f 


i-o 

i<jn 


exp(o‘^j)(m-i)T ) 


(5.5) 


The structured receiver corresponding to (5-5) is denoted 
as structured, receiver A1 and is shown in Pig. 5.2. 
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("b) llYpo t h eeis. t-3stln^ problem A with noncoherent 

l S?P ., 9 . M . ££ i; ♦ In this case the decision rule 

is stated as 


M~1 

I 

m=o 




(5.6) 


This decision rule is simpler to implement than (5 .5). The 
structured receiver corresponding to (5.6) is denoted as 
structured receiver A2 and is shown in Pig. 5*3. 


( c ) Hypothesis testing prohlem B with coherent detection . 
This hypothesis testing problem is solved using the Neyman- 
Pearson criterion for optimisation, A simila.r hypothesis 
testing problem has been solvdXin Section 4*2. Using the 
same procedure, the decision rule can be stated as 


where 


M-l 




X, 


ra»HAX 


m 


MAX 


®m . I 


1=0 



\ exp( jti^(m-i)Tp) 


m > MAX 


(5.7) 


The structured receiver corresponding to (5*7) is denoted 
as structured receiver B1 and is shown- in Fig. 5.4* 


(d) Hypothesis testing problem B with noncoherent 
detection. In this case the decision rule is stated as 



(5.8) 
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Fig. 5. 2 Structured Receiver A1 (Complex Operations). 


r(li 




MAICHED RIMER 
h(t)« f*(TQ~t)^ 
exp(-d ‘^3(^0**^)) 




Sample at 

t = TQ + mTp ; 

m=0,l, . . ,M“1 




MEM WHIOSIRIM 
RIIEER 

MAX ^ ri 

H{f>= E ^ix 
1=0 

exp (~ 3 a ) iTp );-. j 




X, 




**.^1 ■ I 



Rig. 5.3 Structured Receiver A2 (Complex Operations). 
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Close at 
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Pig, 5.4 Structured Receiver Bl( Complex Operations) 


MATCHED PILTBR 

. V .V 


MEM WHITE MI MG 
PILTBR 


h(t)=f (T.-t) X Sample at ^ ~ 

oxp(-jJvt)) 


jllllBnill 

m 


9 


Pig. 5.5 Structured Receiver B2 (Complex Operations) 
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Again this decision rule is simpler to implement than (5*7)» 
The structured receiver corresponding to (5.8) is denoted 
as structured receiver B2 and is shown in Pig. 5.5. 


5 • 2 Performance of the Structured Receiver with Coherent 
D etection 

In calculating the performance of the structured 
receiver I the whitening filter is taken to te just a linear 
filter with transfer function given by (5.2). 


Prom (5.5) and (5.7)> the decision rule is given by 



M-1 ^ ^ % 

y X s*l^ < Dm structured receiver A1 

nio “ “ Hq 


MAX 

where ' exp( j Wjj(m-i ) T^) 

i<.m 


and 



M-1 

I^Iax 




for structured receiver B1 


MAX 

where ^i a)j3(“i*-i)^p) 


m > MAX 


I is a sufficient statistic and is a ccmplex Gaussian 
random variable. In such a case from Section 4.3, we have 


P 


1+^ B[ |I| V%3 - E[ 


(5.9) 
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Define 



(5.10) 

r ie aB defined in (4.13). 

from (5*3) we observe that 

£ » P r" (5.11a) 

where 



(5.11b) 
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E[x x*^) a E[i’ r r*^ 1*^3 = E E[r 
From (5*4) we again observe that 


(5.12) 


"" ^TR ^ 


where 




a 


'MAX 


o. 


\ 

X \ 

\ 


0 


\ 


X 


X 


X ^ 

®M.4X 


\ 


a. 


®r^iE 4:13 = 


DR 


(M-MAZ)xM 

(5.13b) 

(5.14) 


Proceeding along the same lines as in Se-ction 4.3 
and using the same notation for the target and/'interference 
covariance mtrix, we have 

i 

For structured receiver Al» 

F[ |fl VHq] = ^ ® F 4 ® (5.15) 

. • ' 

and 

s[ |i| = 20^ s* P Aj f*® S + 2P“ s** f \ P*'' S ■ 

( 5 . 16 ) 


and 



106 


'a 1 2 j.T - 




For structured receiver Bl, 


( 2h^E^.+NQ)e*^ F F*^ s 

(5.17) 


?b 2 e*^ V A ~ 


CE[ |M7H^3 = 

and 

B[1 M /Hj^] = 20^ 


(5.18) 


+ 2p^ Pjjj Aj fJ^ 


and 


A _ ^t ^IR ^TR ^TR ^!rR _ ^t ^a?R ^IR ^R ^TR 

2P^ ®TR ^TR ^R ^TR ( 2h®E^+NQ) ^"tr'\^Jr®TR 

(5.20) 

Again 2o'^^ and 2p ^ are the expected values of the target 
and interference power in each sample respectively. 




(5.19) 


5 • 3 A Generalised Procedure to Calculate and Pj, for 
Receive re where the Sufficient Statistic is obtained 
hy a Summation of Modulus Squared Correlated Jointly 
Oomnlex Gaussian Random Variables . 

In this section, the expressions for Pj, and Pp for 
receivers with the general configuration shown in Fig. 5.6 
aare derived, ^'(t) is a complex Gaussian rand om process 
under both hypotheses % and respectively and y is a 
sufficient statistic. 
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5(t) 

. - 5 >- 



Pig,5.6 General Configuration of the Receiver 

(Complex Operations) for which Pp and Pj, 
Expressions are obtained. 

We follow the same notation as in Section 5*2, and 
proceed with the derivation using the statistics of 
Under hypothesis let 

^m = = 0,1,...,M-1 ( 5 . 21 ) 

and under hypothesis let 

® ^ » m = 0,1,...,M-1 

( 5 * 22 ) 

where f^'s are jointly complex Gaussian random variables and 
constitute the correlated interference samples and B is a 
complex Gaussian random variable and yE^ B constitutes 
the signal samples. B and I^'s are independent. Therefore, 
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under both hypotheses and are jointly complex 
Gaus&ian random variables. Again following, the same notation 
as in Section 5 . 2 ^ let 


I = [Iq 
and 


~nT~r~ “-t 

■^M-l-' » ^ ~ %-l^ ' [^o ^M“l^ 


f 


B[r] = 0 , E[B] = 0, B[ PI*] = 20^. E[I f h = 2 P" 

(5.25) 

Equations (5»2l) and (5»22) can be written as 


r a VE.|. B u + I 

*** av 

r = I 



where E[(rr*^)/HQ3 = 2P * Aj 


(5.24) 


and E[(r r*^)/H^3 = 2 e| E.j. u u*^ + 2p ® Aj- == 2 y®Aj^ 


2p ^ represents the expected value of the interference power 
in each sample, 2 y^ represents the expected value of the 
combined signal and interference power in each sample when 
the signal is present. 


We now have the following joint density functions, 


and 


_.i expC- r)J 

(2TtP^) Det(Aj) ( 5 . 25 ) 


( 2 ^cy®)^ Itet(Ag) 


exp[. If 


2y' 


(5.26) 
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When no output samples are ignored, the discrete linear 
filter (Whitening, etc.) output may be described by 
means of an MxM transfer matrix 1 such tiiat 

x=fr where x = [\, (5.27) 

When the first M~p output samples are ignored, the required 
filter output may be described by means of a pxM transfer 
matrix such that 




^M-1^ ( 5 • 28 ) 


Por structured receivers A2 and B2 of Section 5.1, P and 
^TR defined in (5.11) and (5*13) respectively. When 

no filter is used, P and Pjij^ are identity matrices. In 
the case of complex weighting of the samples, P and 
aro diagonal matrices with different weights on the diagonal. 


Without loss of generality, we proceed to calculate, 
Pjj and Pp for the case represented by (5.27). 

Then, the joint density function of x ‘s is given by 




— exp[ ^-(x*^ x)3 

( 271^ ^)^ Det(Pj ) 


(5.29) 


and 


^x/I, 


( 2itY*)^ Det(Djj) 


exp[- -^(x*®:^^x) 


2y‘ 


(5.50) 
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where 

Dj = F Aj P ^ and ^ 5*^ (5.31) 

The O'u.'tpu.'t of ijhe &'uimiier/m.'tQgr3.'tor is given by 
M-1 

^ I l^ml“ = ^ (5.32) 

m= 0 

Hermitian laatrices and both are positive 
definite. They may be diagonalised by solving the eigen- 
value problems 

^ ~ ^ ’ k = 0»1,»,,,M— 1 (5*33) 

and 

^R ^R^^ = Pk ^R^^ 5 k == 0,1,...,M-1 (5.34) 

The eigenvalues ^j^'s and Pj^’s will be real and positive, 

and are their corresponding eigenvectors. In 

(5.33) » if all the eigenvalues are distinct, the eigenvectors 
are orthogoml and can be made orthonorml. If there are 
repeated eigenvalues, the corresponding eigenvectors will 
not be orthogonal, but as Dj is Hermitian, these eigenvectors 
will be linearly independent. These eigenvectors can then 

be transformed into orthonormal vectors using the S-ram- Schmidt 
orthonormalisation procedure. Therefore, in general we assume 
that L^^^’s are orthonormal for k = 0,1,.,., M-1. ’'The same 
argument holds true for the eigen values and eigenvectors of 

(5.34) . 



Ill 



Then from elementary matrix theory, 



( 5 . 55 ) 


( 5 . 56 ) 



112 


\fe make the following tranaf ormatione in ( 5 . 32 )! 

^ ^ when Hq is true (5.37) 

and 

Z 7*T ~ 

2 » Lj, X when is true ( 5 , 38 ) 

where z « [z^, z^, 


The suiQraer output in the both cases can be represented by 


y = 



X 



M-1 

J 

ItaO 



(5.59) 


However, the joint density fimction of z when is true 
(P£‘/H^( 2 /Ho) ) is given by 




(2i^P^) 


M 


M-1 

n 

l&xO 


h. 


T M-1 
exp(- -i- y 
2P ^ fc=o 


I 2 


‘ki 


k 


■) ( 5 . 40) 


and the j oint density fmotlon of z wiien Hj is true 



T M-1 

exp( j-J 

2y k=o 



) 


(5. 41) 


The complex Gaussian vector z consists of M independent 
complex Gaussian random variables Zjj-'e. Vhen is true, 
the variance of each random variable is given by 2p ^ 
let 



(5.42) 
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When is truo> will have an exponential density 
function (Prij^/H^C^kAlo)) given by 




Also 


n. 


e.p(. , 0 


2r 


i < 0 


(5.45) 


E 


[exp(-s nj^)/HQ3 * J [exp(-B njj.))p^ 


(1 + 2P Xv s ) 


(5.44) 


It may be observed from (5-59) that 


^y/H. ^y/%) 


0 for y < 0 


(5.45) 


where (y/H^) and Pyyjj^(y/%) are the density functions 

of y given is true and given is true respectively. 

As y/HQ(y given Hq is true ) is obtained by a summation 
of M independent random variables njj.’s, the Laplace transform 
(one sided) of Py^jj (y/HQ) is the product of the Laplace 
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transform of the density functions of Or 




00 


= J«xp(-sy) Py/H^(y/H„)ay = 

(5.46) 


Using (5.44) 

■^y/H 


M-1 

n 


ks=o (1 + 2p^ X. s) 


(5.47) 


Similarly wo can show that 




l&=o (1 + 2 y 3) 


(5.48) 


We then have 


»y/H (y/'io) = 


1^ “ Qxpfsv) ds. 

I M-1 


n (i+2p^ X, s) 

k=o ^ 


Qxp(-S- y) d(s/2P®) 

28 

j . Tpr**= ^ 


-1- 

2nJ 


c-o 


n (1 + s) 

k^o 


and 


1 0 xp(sy) ds 


2110 


J. 


c-o«. M-1 

n (l-+-2Y^ [Xjj. s) 
le=o 


(5,,. 49) 


c+joo exp(-®“ y) d(B/2Y*) 


-Jb— f 
2^^ o-j 


Ml 


j (1 + B) 

kso 


(5.50) 
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The expreeeions to obta,iii Py/u^Cy/H^^) and Py/H^Cy/H^) in 

(5.49) and (5.50) respectively are similar. To simplify 
(5.49) or (5*50), we have to consider the following three 
case s; 


(a) All eigenvalues are distinct. 

(h) All eigenvalues are identical. 

(o) A combination of distinct and repeated eigenvalues. 

Cage ( a ) : All eigenvalues are distinct. 

Wo can express 


1 

M-l 

n (1 + s) 
k=o 


•'^0 , ^ . '^-1 

r ^ # • • ' — 'll ' 'I ,T-rrrr- 

(1 + XqS) (1-+-X;j^s) 1®^ 

(5.51) 


and obtain Aq ... by the partial fraction method etc. 


Then simplifying and taking the inverse Laplace transform 
in (5.49)»we obtain 




M-1 M~1 

ap [ II (1 
itt=o kiso 
k^^m 


Xj^ -1 e3cp(-y/2p 
2P " Xj, 

for y > 0 (5*52) 


Similarly we can obtain 


M-1 M-1 

' Jo f Jo'" 

k^/m 


lifc,,-! e3tp(-y/2Y’ Pn) 

P*)5 


m 


2V‘ |ia 


for y i 0 


(5.53) 



The probability of false alarm and the probability of 
detection are given by 
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F 


and 


M -1 

I * I C n (1 - ^)) exp(T)m/2p 2 X ) 

Mm 0 iiip;0 k=o *ni ' 1' m' 

(5,54) 




<» M-l M -1 n _1 

- I = J - ^)3 eKp(-.,/,ya 

(5.55) 


Case (b); All eigenvalues are identical. 


Then 



Xj^ ss X for k 

= 0,1,. .,,M- 

1 


(5*56) 

and 







« p for k 

=* 0,1,...,M- 

1 


(5.57) 

Now 

from (5*49) we have, 





^y/H 

T 0 + j 00 -1 

dy/H„) = ^ 

f (s+-i- 

-d-§ 

for 

y^ 0 






(5.58) 

Simplifying we have , 






1 

exp(*-y/2p^ 




^y/H, 

(y/Hjj) - jj 

(SP‘ >‘) 

(M~l) 1 


for 

y >0 

(5.59) 

Now 







CO 

1 ‘T y 

exp( 

~y/2p 

^ ) 

F sa 

rp - 

f Pv/ti (y/Ho)dy « 

% ° (2P* 

M 1 

X) ’^T 

(M-1) i 

— dy 



or 
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Pp = 1 - ! 


0 (2pn ) 


2, 


y^*"^ exp(-y/ 2 p = ) 

(M-l) ! 


dy 


(5.60) 


Substituting t = y/ 2 p " and simplifying we get. 


= 1 - r 
0 


.M-1 


(M-l) i 


(5.61) 


uy^(i4-s) . . g 

Now tho integral J jy . ..dy ig Pearson' e form 

o s! 

of the incomplete gamma function and has been defined by 
Pearson [ 58 ] as, 

. »r/ -» \ 


0 bS 


dv 


(5.62) 


Pearson [583hafl also tabulated it upto the seventh decimal 
place for s_< 50. Therefore, frcm (5.61) and (5*62) 

X 


P 


1 - I(- 




, M-1) 


(5.63) 


Following the same procedure, we can obtain 

p , M-1) (5.64) 

D fM 

Case (c) ; A Combination of distinct and repeated eigenvalues. 

M-1 M-1 

In this case l/[ n (1+ V)) or l/[ |T (1+iii, s)) is to 

k=o fe=o 

be expressed as a sum of partial fractions such that the 
methods of Case (a) and Case (b) can be applied to obtain 
Pjj and Pj,, 
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Using tihis procoduTG^ Pj^ and Pp for structursd 
rac^ivors A2 and B2 may be evaltxated. This procedure can 
also bo to evalmte and Pp for receivers using envelope 
square detection followed by post-detection integration. 

Here it must be mentioned that Marcijm [ j) and Swerling [ 8) 
have derived formulae to calculate and P^, for such 
receivers when they are used to detect nonfluctuating and 
fluctuating (Swerling 1-4) targets in white noise, 

Kantor [ 59) has generalised Swerling *s procedure to include 
the effects of clutter for Swerling 1 and Swerling 3 targets. 
Kantor has assumed tho power spectral density of clutter to be 
symmetric. The procedure developed here males no such 
assumptions about the clutter or interference spectrum and 
can be used in the case of both Swerling 1 and Swerling 2 
targets. Also this procedure is easier to use than Kanter's 
procedure which calls for calculating the roots of an Mth 
degree polynomial, M being the number of hits/scan. 

5*4 Performance Evaluation of the Structured Receiver by 
Simulation and Numerical Oomnutation 

The performance of the structured receivers AL, A2, 

B1 and B2 is evaluated by means of computer simulation and 
numerical computation. The performance of the structured 
receiver in all these cases is compared with the performance 
of a moving target indicator (MTI) uslng a double delay line 
canceller after the matched filter. This MTI configuration 
is shown in Pig. 5. 7. 
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Pig, 5.7 MTI Using a Double Iiine Canceller 
( G omple X Ope ra ti one ) . 

Wo have considerod , it appropriate to compare the 
performance of the structured receiver with that of an MTI 
receiver for two reasons. First, the delay line canceller 
MTI receiver is very widely xxsed in operational radar systems 
to distinguish moving targets in clutter. The MTI receiver 
configuration shown in Fig. 5*7 is generally not used and 
the received signal is directly sampled and sent through the 
canceller loops, though radars using pulse compression 
techniq_ues do have some sort of a matched filter front end. 
However, we use a matched filter front end for the MTI 
receiver (whose performance is evaluated here) in conformity 
with the structured receivers of Section 5.1. The second 
reason is that the structured receivers A2 and B2 turn out 
to be similar to an adaptive MTI receiver where the tap 
coofficients (Si's.) and the number of canceller loops (MAX) 
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8,T& s s'fcini&'tGd Edaptivsly from 13116 r6c6iv6d interferenco 
sa^mples* As "the tap coefficients are estimated on the basis 
of combined thermal noiso and clutter* it is conjectured that 
the performance of the structured receiver will be superior 
to that of an adaptive MTI receiver where the Cffsjxceller loops 
are designed to combat clutter alone. 

The perfoimance of the MTI receiver is evaluated for 
the following cases using the procedure discussed in 
Section 5.5. * 

(a) When all the M pulses at the output of the canceller 
are taken into account. This receiver is denoted as 
MTI receiver A. 

(b) When the initial two pulses at the output of the 
canceller are ignored. Here the Bvunmation in Mg. 5.6 ,, 
is from m = 3 to M-1. This receiver is denoted as 

MTI receiver B. • 

As in Chapter 4, the performance evaluation is done 
for varying conditions of target power* target velocity, 
clutter power, clutter spectral width, thermal noise power 
etc. Again, the clutter returns are assumed to have a 
Gaussian shaped power spectral density. 

The entire simulation scheme is discussed in 
Appendix A* As in Chapter 4, zero-mean complex Gaussian 
clutter samples with Gaussian shaped power spectral density, 
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aro generated on the computer. Zero-mean complex Gaussian 
white noise samples are also generated on the computer and 
these are added to the clutter samples to obtain interfeiance 
samples These Ij^'s correspond to the output of the 

ma.tched filter in Figs. 5.2 - 5.7 when no target is present. 

Then, the complex foim of the maximum entropy method is used 
to determine the tap coefficients (\’b) and order (MAX) of 
the whitening filter. The performance evaluation of the 
structured receivers Al, A2, B1 and B2 and MTI receivers 
A nnd B ie dons using the expressions and formulae 
obtained in Sections 5.2 and 5 •3* This evaluation is 
carried out for a large number of cases by varying Pg,, 

Nq, fp and f^ where these terms are as denoted in Chapter 4-. 

The receiver operating characteristics (ROC) for these cases 
are given in Figs. 5.8 “ 5*20. 

The following observations on the overall performance 
of 'the structured receivers Al, A2, B1 and 32 f and on their 
comparative performance' versus the MTI receivers A and B, are 
made : 

(a) In the situa-tions considered, the performance 
of the structured receiversCwitheitler coherent or noncoherent 
detection) based on hypothesis testing problem B is far superior 
to that of MTI receiver B. likewise, the performance of the 
structxrred receivers (again with eithei'' -coheront: or noncohoi^mt 
detection) based on hypothesis testing problem A’ is superior 
or at least equal to that of MTI receiver A. 
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(Td) When there is no clutter, the performance of the 
structured receivers tesed on hypothesis testing problem A is 
slightly better th'^n the performance of the structured receivers 
based on hypothesis testing problem B« However, the performance 
of the structured receivers based on hypothesis testing problem 
B is far superior to that of the structured receivers based on 
hypothesis testing problem A even when the clutter power (B ) 
ia only 20 dB above the white noise power (N^) and 10 dB above 
the target power (P^). This superiority in performance 
increases with an increase in clutter power. Likewise, under 
similar clutter conditionsj the MTI receiver B performs better 
than the MTI receiver A, At a target Doppler shift of 100 Hm 
and Pj equal to 10, the performance of the structured receiver 
A2 and the MTI receiver A approach the Pj^ci Pj, limit when the 
clutter power is 50 dB above the white noise power. 

(c) As ey.pected, the structured receivers with coherent 
detection perform better than the structured receivers with 
noncoherent detection, 

(d) In general, the performance of the structured receivers 
improve as the clutter spectral vridth decreases. 

(e) The performance of both the structured and the MTI 
receivers is dependent on the target Doppler shift. For 

P « 1000 and f, « 1 Hz, the performance of the MTI receiver 
B deterio»t©s sharply when the target Doppler shift decreases 
from 100 Hz to 40 Hz. However, the performance deterioration 



is much less in. the case of the structured receivers B1 and 
B2 for the same change in target Doppler shift* In the case 
of the structured receivers B1 and B2» this performance 
deterioration, with a decrease in target Doppler shift, 
is more rapid as the clutter spectral width increases* 

(f) It ie observed that when f^ is equal to 10 Hz, 
the performance of the MTI receiver B is thermal noise 
limited upto clutter power levels of 40 dB above the 
thermal noise power. This is obviously not so in the case 
of the structured receivers. 

The performance of these structured receivers is 
inferior to the performance of the structured optimum 
receivers of Chapter 4* However, the implementation of the 
structured receivers (particularly with noncoherent detection) 
is easier than tho implementation of the etructuiod optimum 
receivers. 
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CHAPTER 6 

ooNoiupioyrs! 


In tihis thssisj an. opljimuni radar recaiver for 
detecting a Swerling 1 target in the presence of clutter 
is dorilwd' and it is found that implementing this re ceiver 
in a changing clutter environment is extremely difficult 
with present day technology. In view of this, this thesis 
primarily deals with the design and study of two adaptive 
sub optimum receivers for detecting a Swerling 1 target in 
the presence of clutter. These receivers are adaptive in 
the sense that they estimate the required statistics of 
clutter and adjust or update their parameters in response 
to ctongOB in the statistics of clutter. This is achieved 
through the maximum entropy method (MfiM) of spectral ana.lysiB 
of the received signal during those scans when there is no 
target or a decision to that effect is made. 

In both the suboptimum receivers for detecting 
a Swerling 1 target, a matched filter (matched to a single 
pule© return) structure is imposed. A matched filter 
structured receiver is relatively easy to implenient and 
is indeed the optimum receiver in the absence of clutter. 
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The first subojptimm receiver, termed Structured 
Optimum Receiver, is an optimum receiver within the constraints 
by its matched filter structure, A study is made of 
how tht- MiM could be used to estimate the required clutter 
EtdistiCB during those scans when no target is present in 
the rangtJ bin of interest, thereby enabling the para.nieters of 
the structured optimum receiver to be determined adaptively, 

Tlw performance of such a receiver is evaluated for both the 
cases vdien the required clutter statistics area priori known 
and when, the MfilM is used to estimate the required statistics 
of clutter. Performance evaluation studies indicate that the 
Mi:Hj is able to estiraate the required clutter statistics 
'reasonably well' with only 16 samples. In general, the 
performance of the structured optimum receiver improves as 
the clutter spectral width decreases and the target Doppler 
shift increases, At a target Doppler shift of 100 Hz or 
more, a pulse repetition frequency (PEP) of 540 Hz and a 
clutter spectral width of 10 Hz or less, the performance 
of tlio structured optimum receiver approaches the maximum 
attaiimble limit (when no clutter is present) to within 
about 20 porcont (in terms of signal to interference ratio) 
oven when the clutter power is as high as 50 dB above the 
thormnl noise power. In general, the performance of the 
structured optimum receiver is far superior to that of a 
conventional matched filter receiver. However, at low target 
Doppler shifts (particularly at 20 Hz. or less), the performance 

of both the receivers is poor. 
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su'boptiinum receiver is easier to impleinent than the first 
sub optimum receiver. 

A study of the performance of both the structured 
optim\M and structured receivers carried out at other PRI''s - 
though pot reported here indicates trends similar to those 
presented in the thesis. 

In the procedure followed to estimate the required 
clutter statistics by MEM, an averaging (as described in 
Appendix A) over twenty five scans (when no target is present) 
is done. This procedure ignores the errors due to missing a 
target during a particular scan and considering the returned 
signal in that scan to estimate the clutter statistics. To 
take this error into account, one has to consider the a priori 
probability of a target being present at the given range bin 
during any particular scan, as well as the probability of 
miss of a target. In practice, this apriori probability of 
a target being present is extremely low, and therefore it 
is felt that the errors due to missing a target will not be 
significant when an averaging over a sufficient number of 
scans is donei 

The performance figures of the structured optimum 
receiver indicate that its performance is 'reasonably good' 
except at low target Doppler shifts and large clutter spectral 
width#. At low target Doppler shifts, the performance is 
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©xtremely poor* To obtain better perforioance at low 
Doppler shifts ^ it is necessary to search for other 
implementable receiver structures. One such receiver 
structure (in the notation. of the earlier chapters) that 
could be tried is shown below. 


^tt) 


at 


M ~1 ^ * 

Jlq Gxp(-3aij5t) 



M-1 1 ^ 

^ mlo hlo”^^ ®^P(-3“D“^p) f*(t-kT-inTp) 


This has tho same structure as the optimum receiver except 
that tho coefficients in the optimum receiver have been 

replaced by exp(jw^mTp) terms, could be chosen such 

that tho clutter component in the received signal is completely 

cancelled. Using the notation of the earlier chapters, 

L ^ 

T ^mk f(t-kT-mT„) is the clutter signal due to the mth 
pulse, f(t-mT ), in the transmitted pulse train; When this 
return is correlated with g*(t), the output due to the clutter 
component alone can be written as 


Ij I* 1 

V oxp(-ja,i,-Tp) V 






where f(t-kT)f*(t) exp(-3WDt) dt 
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To make this clutter compoviant equal to zero, we can put 
the above expression to zero and solve for Or 

iC 

■ ^ - 

I i((k-n)T,0) = ~(p(kT, - 0)^) 

for k = 0,1,... ,1 

We observe that these Pj^*s do not depend on m and choosing 
Pjj.’s in such a fashion would result in the complete 
cancellation of clutter. However, the power due to the 
white noise component in the received signal increases at 
the output of the receiver. One way to compensate for this 
is to increase the transmitter power or the number of hits/ 
scan. Another way to choose Pjj.*s is to maximise the signal 
to interference ratio and perhaps use a gradient algorithm 
to recursively compute Pjj.'b. 

The performance evaluation of the various receivers 
ie done by considering the statistics of the samples at 
the output of the matched filter. Ho consideration was 
given to the design of the basic signal waveform (f(t)) 
to reduce the effects of clutter. By a proper choice of 
the basic signal waveform (f{t)), the clutter power (2h^) 
at the output of the matched filter itself could be reduced 
and this would inprove the performance of the structured 
optimum and structured receivers considerably. Thus if 
f(t) is BO chosen that its ambiguity function (©(t, w)) has 
a low value at small values of ^ , the performance of these 
receivers at low Doppler shifts would improve considerably. 
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PBRPORMAIQE EVALUATIOM' BY COMPUdER SIMUIATION AM) 
MMERICAL aOMPUTAHONS 


The performance of the structured optimum, structured, 
matched filter and moving target indicator (MTI) receivers 
were evaluated for Swerling 1 target models by computer 
simulation and numerical computations. The computer used 
was the DiilC-1090 system. This appendix describes the 
simulation scheme, A flow cha,rt' describing the 
recursive algoritiim for the complex form of the maximum 
entropy method (MEM) is also included. 

The objective of the computer simulation was to 
evaluate the effectiveness of the MEM in estimating the 
required total interference statistics in the case of the 
structured optimum receiver»and to estimate the required 
whitening filter in the case of the .structured receiver. As 
is commonly done in most radar performance evaluations, the 
clutter returns were assumed to have a (jaussian eha.ped power 
spectral density. The number of hits per scan (n) of the 
radar was taken to be 16. 

Complex Gaussian clutter samples were generated 

on the computer using the method described in [ 603 . Complex 
white Gaussian noise samples (w^^’s) were also generated on the 
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computor and tiiese were added to to obtain the total 

interference samples (ij^'s). These interference samples 
correspond to the matched filter ontpat for a particular 
range bin during those scans when no target is present. 

Then the complex form of the M&iM was used to estimate the 
inverse covariance matrix (Aj^) of the interference samples 
in the case of the structured optimum receiver and to 
estimate the whitening filter tap coefficients {^■^) aud its 
order (MAX) in the case of the structured receiver. ‘As 
indicated in Chapter 2> the inverse covariance matrix (A^ ) 
can "be estimated without the need to estimate Aj and then 
invert it. Also for the reason stated in Chapter 2# the 
maximum order of the autoregressive (AR) model being fitted, 
on the samples was restricted to 8, 

——I 

In order to improve the estimates of Aj- and the tap 
coefficients of the whitening filter, twenty five different 
realisations of the complex Gaussian interference samples 
(corresponding to those scans when no target is present) 
were generated on the computer and the following averaging 
scheme used: 

(a) Th'Om each realisation of 16 samples, the MEM was 
used to calculate the AR model coefficients for . 

m « 0,1,,..,8 and k - 0,1,... ,« by solving equation (2.44) 
reproduced below in terms of complex envelopes, 
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m ^ „ 

%k 2 ~ 0,1,. . . ,iti 

®^5E 0ylj.**>8 

4W 

= 1 for aiDy ja. 

m denotes the order of the ilR model being fitted, R(k) ' s 
de^noto the element of the interference covariance matrix 
( Ai ) and P^ denotes the prediction error power for a 
prediction error filter of length m+l. These equations are 
solved using Burg’s procedure outlined in Chapter 2* The 
final prediction error (;pPE(m)) as stated in (2.55) was 
also calculated for 1,2,...,8, A flow chart 
of the MEM recursive algorithm for a complex -valued time 
Berios is given in Pig. Al.l. The notation used in this 
figure is that of Chapter 2, only the subscripts are 
enclosed in parenthesis. MAXOED in Pig. Al.l indicates the 
maximum order of the AR model for which the search is to be 
carried out. In our case MAXOED is set to 8, 

(b) These and PPI(m) wrere averaged over twenty 

five different realisations of the interference samples. Then 
AkailcO’s minimum PPB criterion was used to select the /iR model 
order (m s MAX). Thon was calculated using. (2-59). The 
whitening filter coefficients are given by “jj-MAX 
k s 0,1,. . , ,MilX. 

The performance of the structured optimum, 
structured, MTI and matched filter receivers was then 
evaluated by using the formulae and procedures 
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described in Chapters 4 and 5. The performance of these 
roc^-ivers was evaluated for a large number of cases by- 
varying the clutter power, clutter spectral width, target ' 
power, target velocity and thermal noise power. Appendix ® 
gives a listing of tho main program and the subroutines used 
in the performance evaluation. 
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